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ABSTRACT 


AlGaN  DIRECTIONAL  COUPLER  SWITCH 

Crespo  Antonio 
University  of  Dayton,  2002 

Advisor:  Dr.  Mohammad  A.  Karim 

This  dissertation  describes  an  attempt  to  develop,  including  material  processing,  a  zero-gap 
directional  coupler  using  the  AlGaN  material  system  in  an  n-i-n  configuration  that  was  grown 
by  MBE.  This  work  compiles  the  optical  and  electro-optical  properties  and  models  the 
behavior  of  AlGaN  integrated  optic  waveguides.  The  switches  are  based  on  the  zero-gap 
electro-optic  directional  coupler.  The  models  are  compared  with  the  experimental  results  that 
were  obtained.  The  models  used  were  the  Marcatili  and  the  beam  propagation  method.  Models 
for  predicting  the  index  of  refraction  of  AlGaN  versus  A1  concentration  are  treated.  An 
advantage  of  the  AlGaN  integrated  optic  waveguides  and  switches  is  the  potential  to  work  at 
wavelengths  from  IR  to  the  blue  region.  AlGaN  has  the  potential  to  work  at  high  temperatures 
and  in  harsh  environments.  The  wavelength  chosen  for  this  work  is  632.8  nm  from  a  HeNe 
laser. 

The  results  show  that  an  electro-optic  switch  is  viable  in  the  AlGaN  system.  Material 
maturity  is  necessary  to  improve  performance.  The  models  chosen  were  shown  to  be 


inaccurate  by  themselves.  Together  they  are  a  viable  source  for  the  development  of  a 
directional  coupler  surrounded  by  air. 
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CHAPTER  I 


INTRODUCTION 


As  technology  advances,  electronic  circuits  increase  in  complexity  and  performance. 
This  increase  also  requires  an  increase  in  the  number  of  interconnections  per  chip  and  an 
increase  in  the  interconnection’s  speed,  otherwise  the  bottlenecks  created  will  limit  the 
overall  system  performance.  The  interconnections  include:  chip-to-chip  on  a  common 
multichip  module,  multichip  module  to  multichip  module,  and  board-to-board 
communication  through  a  back  plane.  The  current  trend  in  technology  is  heading  towards 
multi-gigahertz  communication  between  processing  circuits  and  greater  than  2000  pins 
per  chip  in  the  near  future.1  The  increasing  cost  in  terms  of  power,  size,  component  count 
weight  and  reliability  of  implementing  the  required  interconnects  electrically  has 
motivated  serious  consideration  of  replacing  electrical  interconnects  with  optical 
interconnects.  The  data  rates  at  which  cost-effective  communication  through  an  electrical 
interconnect  can  be  transmitted  is  limited  to  low  (in  the  10’s  of)  gigahertz  rates.  ’  This 
results  from  the  fact  that  electrical  wires  and  striplines  incur  frequency  dependent 
dissipative  losses  due  to  conductive  and  dielectric  properties.  Electric  lines  have 
frequency  dependent  parasitic  coupling  losses  due  to  radiative  and  surface  wave 
propagation  properties.  However,  optical  interconnects  can  transmit  data  at  higher 
frequencies  with  orders  of  magnitude  lower  losses  than  their  electrical  counterparts. 
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MCM  substrate 


(a)  (b)  (c) 

Figure  1.1  Packaging  interconnection  levels,  (a)  level  3  (wirebond),  (b)  level  2  (TAB), 

and  (c)  level  1  (C4);  from  reference  [4] 

In  chip-to-chip  connection  bonding,  large  pads  are  needed  to  make  connections  between 
the  chips  and  the  MCM  substrate.  The  large  capacitance  and  inductance  of  the  pads  tends 
to  increase  the  interconnect’s  power  dissipation  and  reduce  the  interconnect’s  bandwidth. 
The  number  of  contacts  needed  to  connect  the  chip  to  a  substrate,  known  as  contact  level, 
also  influences  the  bandwidth,  Figure  1.1  shows  some  examples.  Transmission  line 
effects  will  dominate  the  interconnection  when  line  distances  are  greater  than  1/1 0th  of  the 
electrical  signal  wavelength.  This  results  in  attenuation  of  the  electrical  signals  through 
skin  effects,  degradation  of  signal  integrity  due  to  reflections,  and  increased  power 
dissipation  occurring  with  terminating  resistors.5  Within  the  MCM  the  electrical 
interconnects  are  well  suited  to  handle  communications  below  the  10’s  of  GHz  and  no 
significant  increase  in  MCM  size.  For  the  MCM-MCM  communications  the  bottleneck 
increases.  The  bandwidth  and  the  connection  densities  are  reduced  by  more  than  an  order 
of  magnitude.  In  the  backplane  level,  electrical  interconnects  would  have  difficulty  in 
meeting  a  bandwidth  times  length  times  density  product  exceeding  20  GHz  cm/mm. 

Optical  communication  has  an  inherently  high  bandwidth.  The  optical  carrier  has  a 
frequency  of  the  order  of  ~1014  Hz,  orders  of  magnitude  higher  than  any  modulation 
source  available  today.  Thus,  it  is  able  to  meet  the  communication  bandwidth  for  the 


3 


Figure  1.2  Incidence  at  an  interface 


foreseeable  future.  Photons  will  propagate  across  transparent  material  with  little 
attenuation  and  power  dissipation.  The  attenuation  and  power  dissipation  will  not  vary 
with  the  modulation  frequency.  By  varying  the  index  of  refraction  of  the  transparent 
material,  light  can  be  guided  along  a  predetermined  path  (optical  interconnection)  to  its 
final  destination.  This  is  due  to  the  total  internal  reflection  (TIR)  phenomena.  When  light 
is  incident  at  an  interface,  see  Figure  1.2,  it  is  transmitted  and  reflected  according  to: 

n{  sinG;  -nt  sin0t  Snell’s  law  (  1-1 ) 

0;  =  0r  Reflection  law  (  1-2 ) 

When  nt  <  »,  equation  (  1-1  )  tells  us  that  there  is  a  critical  angle  0 ,  =  0  c  for  which 
0t  =  90°.  At  this  angle  TIR  occurs  and  light  is  totally  reflected  back  into  the  medium. 
Therefore  if  a  higher  index  media  is  sandwiched  between  two  lower  index  media,  light 
can  be  confined  to  the  lower  index  media. 

Current  optical  fibers  have  attenuation  in  the  order  of  0.2  dB/km  and  integrated  optical 
interconnects  are  in  the  order  of  0.2  db/cm,  for  plastic  waveguides.  The  advantages  of 
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optical  interconnects  include  the  difficulty  of  photons  interacting  with  each  other,  their 
immunity  to  electro-magnetic  interference  (EMI),  increased  I/O  density,  increase 
communication  speed,  reduced  crosstalk,  and  their  ease  of  impedance  control.5,6  Optical 
interconnects  are  starting  to  be  utilized  in  chassis-to-chassis  and  board-to-board 

O  a  O 

communications  using  glass  optical  fibers,  although  with  some  limitations.  ’  ’  These 
optical  fibers  are  lightweight  and  need  little  shielding  from  cross  talk.  For  a  rapid 
introduction  of  optical  interconnects  into  the  computer  arena  its  manufacturing 
procedures  must  be  compatible  with  the  electrical  components,  for  which  a  huge 
industrial  base  is  already  in  place.  The  point  at  which  optics  will  become  competitive 
with  electronics  is  a  function  of  cost,  interconnection  distance,  signal  density,  bit  rate,  and 
fanout.4  If,  in  addition,  switching  operations  could  be  incorporated  into  the  optical 
interconnection,  without  optical-electrical-optical  conversion,  reconfigurable 
interconnects  would  be  achieved  without  a  reduction  in  the  available  bandwidth. 
Aluminum  Gallium  Nitride  (AlGaN),  and  other  III-V  materials,  possesses  an  electro- 
optical  effect  that  can  be  utilized  to  fabricate  optical  switches  and  modulators. 

AlxGai_xN  has  recently  enjoyed  an  increased  interest  because  of  its  wide  and  direct  band 
gap,  which  ranges  from  3.4  -  6.2  eV  for  x  ranges  of  0  -  1  for  the  wurtzite  structure.  This 
makes  it  ideally  suited  for  blue  to  UV  lasers.  An  additional  advantage  of  this  material  is 
that  electronic  devices  made  from  it  are  more  tolerant,  than  Si  or  GaAs  devices,  to 
elevated  temperatures  and  caustic  environments.9  Also,  GaN  possesses  a  higher  velocity 
of  saturation  (~2.5  x  10  )  than  GaAs,  enabling  it  to  operate  at  high  speed,  and  a  lower 
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dielectric  constant  (er  «  9)10  than  Si,  which  reduces  its  capacitive  contribution.  Because 
AlxGai.xN  has  not  been  grown  in  bulk  crystals  big  enough  to  form  substrates,11  alternate 
materials  with  lattice  constant  close  to  it  have  been  used  as  substrates.  Only  recently  has  it 
been  possible  to  overcome  some  of  the  lattice  mismatch  problems,  by  using  an  optimized 
buffer  layer  on  sapphire.10  Lattice  mismatch  between  the  substrates  and  the  AlxGai_xN 

9  2  12 

layer  is  a  major  contributor  to  a  high  defect  density  (~5  x  10  defects/cm  ,  for  AI2O3 ) 

O 

with  a  corresponding  optical  loss  of  ~10  dB/cm,  but  researchers  perceive  today  that  10 
defects/cm  should  be  realized  in  the  near  future,  if  the  current  effort  in  the  technology  is 
maintained.  With  more  effort  placed  on  reducing  the  defect  density  and  some  additional 
time,  a  reliable  defect  density  of  10  defects/cm  should  be  obtainable.  A  major 
reduction  in  defect  density  should  be  achieved  when  AlxGai_xN  substrates  are  realizable. 
Additionally,  present  growth  techniques  have  high  n- type  background  carrier 
concentrations.  It  is  difficult  to  dope  GaN  p- type,  but  that  is  beginning  to  be  overcome. 

Electro-optic  modulation  consists  of  an  electro-optically  active  dielectric  medium  and  a 
means  for  applying  an  electric  field  to  it.  Application  of  an  electric  field  causes  a  change 
in  the  refractive  index  of  the  dielectric  medium.14  The  relationship  between  the  electric 
field  and  index  change  for  an  electro-optic  crystal  is  given  by  the  series  tensor  equation:15 

(1-3) 

where  E  is  the  applied  electric  field  vector,  r  is  the  linear  electro-optic  tensor  and  5  is  the 
quadratic  electro-optic  tensor.  Terms  higher  than  quadratic,  in  the  expression  above,  are 
neglected  due  to  their  small  contribution  to  the  index  change.  In  crystals  where  the  linear 


A  — r  =  r,,,E.  +  s-^E.E,  +■■■ 
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electro-optic  effect  is  present  the  quadratic  effect  is  also  neglected.  Using  only  the  linear 
portion  and  taking  advantage  of  crystal  symmetry,  equation  (  1-3  )  reduces  to: 


(1-4) 


in  which  contraction  notation  is  used,  see  reference  15. 

AlxGai_xN  possesses  an  index  of  refraction  range  from  -2.35  -  2.0  (for  x  =  0  -  1)  which 
can  be  tuned  with  alloy  composition.  This  index  range  makes  it  a  close  match  with  that  of 
communication  grade  optical  fibers.  The  advantage  is  reduced  reflection  of  the  optical 
signal  due  to  any  index  mismatch.  Additionally  GaN  possesses  an  electro-optical  effect 
that  can  be  utilized  to  fabricate  optical  switches.  The  measured  electro-optic  coefficient 
are  r33  =  1.91  ±  0.35,  r3J  =  0.57  ±0.11,  and  r5\  =  3.1  ±  0.2  pm/V  at  633  nm,16’ 17  which  are 
larger  than  AlGaAs  .  Using  these  parameters  and  the  index  ellipsoid  for  wurtzite  GaN 
([6mm]  point  group  symmetry)  under  an  applied  electric  field  E  =  [Ex,  Ey,  Ez]  given  by: 17 


2  rn^z 

nn 


2  ^13  Ez 

n„ 


+  —  +  r33  EZ 

n„ 


+  2yzr51E  +  2 xzrslEx  =1  (  1-5  ) 


we  can  determine  the  index  of  refraction  that  an  optical  beam  is  subjected  to.  This 
information  is  necessary  for  optical  switching. 

The  index  change  causes  a  change  in  the  optical  characteristics  of  the  material,  which 
can  be  exploited  to  manipulate  the  propagation  of  light  through  the  material.  With  this, 
we  can  create  different  types  of  integrated  optical  switches  (modulators)  that  can  be 

O  1  O 

fabricated  in  a  single  waveguide  structure  ’  .  Among  them  are: 
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Figure  1.3  Basic  electro-optic  switching  structure 
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1)  Polarization  switch  ’  :  This  switch  works  by  changing  the  polarization  of  the 
incoming  optical  signal,  and  using  a  polarization  sensitive  detector  (analyzer)  to 
determine  the  state  of  the  signal.  The  simple  single  channel  switching  structure  on  Figure 
1.3  can  be  used  to  accomplish  this.  Launching  an  optical  beam  linearly  polarized  at  45°  to 
the  x  and  v-axes  and  applying  a  voltage  will  make  the  polarization  rotate.  This  is 
accomplished  by  the  fact  that  the  phase  changes  only  for  the  waves  polarized  in  the  y 
direction  and  not  in  the  x  direction.  Difficulty  in  fabricating  an  effective  analyzer 
monolithically  has  limited  the  use  of  polarization  switches. 


V 


Figure  1 .4.  Mach-Zehnder  interferometer 
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2)  Phase  switch:  By  applying  a  voltage  to  the  device  in  Figure  1.3  an  index  change  will 
be  induced  in  the  device  for  an  incoming  light  polarized  along  the  y  direction,  thus 
changing  the  phase  of  the  polarized  light.  In  order  to  detect  this  phase  change  a  coherent 
detection  system  must  be  used.  An  integrated  Mach-Zehnder  interferometer,  Figure  1.4, 
is  such  a  system.  Using  the  Mach-Zehnder  with  a  n  phase  shift,  the  incoming  light  can  be 
switched  on  and  off  to  create  a  digital  signal. 

3)  Directional  coupler  switch:18,19  When  two  waveguide  channels  are  closely  spaced 
they  will  function  as  a  directional  coupler,  in  which  the  energy  is  synchronously 
transferred  from  one  guide  to  the  other.  In  order  to  obtain  a  100%  power  transfer  the 
waveguides  need  to  be  identical.  The  directional  coupler  can  be  made  into  a  switch  by 
adding  electrodes,  as  shown  in  Figure  1.5.  The  application  of  a  voltage  will  change  the 
index  of  refraction  of  the  electro-optic  material.  This  will  change  the  propagation  constant 
of  one  of  the  waveguides.  With  an  appropriate  voltage,  the  waveguide  will  switch  guides 
from  which  the  optical  energy  is  coming  out. 


Figure  1 .5  Active  Directional  coupler 
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The  directional  coupler  switch  will  be  the  focus  of  this  dissertation,  specifically  the  case 
were  the  gap  (g  in  Figure  1.5)  between  the  waveguides  is  zero,  zero-gap  directional 
coupler.  The  relationship  between  length  and  voltage  necessary  to  switch  100%  of  the 
light  from  one  waveguide  to  the  other  (both  identical)  is  given  by: 


Xd 


(1-6) 


where  X  is  the  free  space  wavelength,  n  is  the  index  of  refraction  of  the  interaction  region, 
r±  is  the  electro-optic  coefficient,  and  VK  is  the  voltage  required  for  complete  switching. 
Using  this  equation  with  the  following  parameters:  Vn  =  5  v;  X  =  830  nm;  d  =  1  |im; 
rjk  =  i*33  =  1.91;  n  =  2.316  [Ref.  16];  we  get  a  GaN  zero-gap  directional  coupler  switch 
with  L  =  7  mm. 

The  advantage  of  investigating  a  directional  coupler  is  that  it  can  perform  a  number  of 
useful  functions  in  thin  film  devices.  These  functions  include  power  division,  modulation, 
switching,  frequency  selection,  and  polarization  selection.20  Additionally  the  zero-gap 
directional  coupler  uses  less  voltage  to  switch  than  conventional  directional  couplers.3 
Furthermore,  using  AlGaN  instead  of  AlGaAs  to  make  these  devices  will  have  the 
following  advantages: 


Closer  index  of  refraction  match  to  the  optical  fibers  and  polymer  waveguides 
Potential  for  working  with  wavelengths  from  the  IR  to  the  blue  region 


Tolerate  harsh  conditions 
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Directional  couplers  have  been  around  for  some  years,  especially  in  LiNb03,  and  are 
well  understood.  The  concept  of  directional  coupler  as  a  modulator  (switch)  was  proposed 
as  early  as  1969  by  Marcatili,  but  quite  a  few  years  passed  before  a  working  device  was 
realized,  largely  because  of  the  difficulty  of  fabricating  the  dual-channel  structure  to  the 
required  tolerances.  Somekh  fabricated  dual-channel  directional  couplers  in  GaAs  with 
100%  coupling,  and  theoretically  analyzed  the  case  of  non- identical  waveguides  for  a 
directional  coupler  in  1973-1974.  Taylor  theoretically  analyzed  the  performance  of  a 
directional  coupler  with  three  electrodes,  see  Figure  1.6.  Campbell  et  al.  were  the  first 
ones  to  produce  an  operational  directional  coupler  with  three  electrodes  in  GaAs,  in 
1975.  Papuchon  et  al.  reported  the  first  two-electrode  operational  directional  coupler 
modulator  in  Ti:LiNb03,  around  1975.  In  1977  Papuchon  et  al.  introduced  the  concept 
of,  and  demonstrated,  a  directional  coupler  with  no  gap  between  the  waveguides,  a  zero- 
gap  directional  coupler,  in  Ti:LiNbC>3.  More  recently,  Grote  explored  zero-gap 


Figure  1.6  Three  electrode  directional  coupler 
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directional  coupler  E-0  switches  in  AlGaAs.  For  this  dissertation  a  zero-gap  directional 
coupler  (sometimes  called  BOA )  in  AlGaN  will  be  attempted. 

The  objective  of  this  research  was  to  investigate  the  behavior  of  an  optical  directional 
coupler  switch  fabricated  from  AlGaN  in  an  n-i-n  configuration.  The  research  explored 
the  possibility  of  the  AlGaN  material  system  as  the  basis  for  integrated  optical  devices, 
other  than  laser  diodes.  AlGaN  is  capable  of  withstanding  elevated  temperatures  and 
harsh  environments  and  has  the  potential  of  harboring  high  speed  high  power  integrated 
circuits.9,10  This  research  establishes  what  areas  of  improvements  are  required  of  the 
material  system  for  integrated  optical  devices.  In  addition,  these  optical  devices  could  be 
designed  to  operate  in  the  region  from  UV  to  IR.  This  investigation  also  developed  the 
necessary  processes  and  techniques  required  to  create  electro-optical  structures  in  this 
new  material  system.  AlGaN  characteristics  and  requirements  are  very  different  from  the 
current  material  systems  used  in  the  industry  ( i.e.  AlGaAs,  Si,  LiNbCh,  etc.  ).  These 
structures  have  features  that  are  both  very  long  and  thin,  which  contrasts  with  electronics 
where  both  dimensions  are  comparable  in  size.  The  research  also  determines  the 
refraction  indices  from  published  data  of  the  AlGaN  material  in  a  layered,  electro- 
optically  active,  device.  The  n-i-n  configuration  was  used  to  concentrated  the  electric 
field  across  the  core  region 

The  directional  coupler  can  be  used  as  a  switch  for  reconfigurable  optical  interconnects 
and  digital  circuits,  such  as  network  routing  switches,  without  the  optical-electrical- 
optical  conversion.  This  is  applicable  to  parallel  systems  with  high  speed  interconnects 
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which  could  be  reconfigured  by  software  depending  on  the  application.  The  use  of  the 
zero-gap  directional  coupler  will  enable  the  use  of  a  simple  configuration  for  obtaining 
complete  power  transfer,  theoretically,  with  a  single  bias  voltage  per  device.  Thus 
reducing  the  complexity  of  the  circuit.  The  experimental  directional  coupler  will  be 
compared  with  the  results  from  the  Marcatili  method  and  the  beam  propagation  method  to 
determine  their  suitability  to  predict  the  behavior  of  the  structure. 

Dissertation  Organization 

An  overview  of  the  project  and  reasoning  has  been  presented  in  this  chapter.  The 
following  chapters  will  elaborate  further  into  the  topics  and  recount  in  chronological 
sequence  the  experimentation.  Chapter  II  will  review  the  theory  of  the  optical  waveguides 
as  well  as  the  Marcatili  and  the  beam  propagation  methods.  The  electro-optic  effect  and 
the  non-linearity’s  of  crystals  will  be  reviewed  in  Chapter  III.  The  directional  coupler 
theory  will  be  reviewed  in  Chapter  IV.  The  current  status  of  the  AlGaN  material  system  is 
described  in  Chapter  V.  A  chronological  exposition  and  description  of  the  experiments 
and  results  is  found  in  Chapter  VI.  And,  in  Chapter  VII  a  discussion  of  the  results  and 


conclusions  will  be  made. 


CHAPTER  II 


WAVEGUIDES 


As  mentioned  in  the  previous  chapter  light  can  be  confined  using  TIR.  This  is 
accomplished  by  surrounding  a  material  that  has  a  high  index  of  refraction  with  a  lower 
one.  Taking  the  simple  case  of  a  slab  waveguide,  see  Figure  2.1,  light  is  launched  at  the 
interface  between  air  and  the  wave  guiding  core.  For  light  to  be  confined  in  the  core 
region  the  angle  that  the  light  ray  makes  with  the  normal  to  the  core-cladding  interface 
must  be,  from  Equation  (  1-1  ). 


(2-1) 


That  would  ensure  that,  at  the  core-cladding  interface,  there  would  be  TIR.  Otherwise, 
light  would  be  escaping  through  the  cladding  region.  Related  to  this  angle,  there  is  an 


external  angle,  0a,  which  defines  the  radius  of  the  cone  of  possible  incident  angles  at  the 


air-core  interface  for  TIR.  This  angle  is  used  to  determine  the  numerical  aperture  (NA)  of 
the  fiber  with  the  following  equation: 


Figure  2.1  Slab  Waveguide 
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NA  =  na  sin  0a  (  2-2  ) 

That  is  widely  used  to  describe  optical  fibers.  Not  all  of  the  rays  incident  at  angles  below 
0a  will  propagate  through  a  slab,  only  those  whose  phase  fronts  coincide.  This  could  be 
more  easily  seen  in  Figure  2.2.  The  dashed  lines  represent  the  ray’s  phase  fronts.  The 

phase  fronts  will  coincide  when  the  total  phase  change  incurred  by  the  segment  CD  is  a 

22 

multiple  of  2%.  This  occurs  when  the  following  relation  is  maintained: 

nx  (s2  —  Sj )  k0  +  <f>2  +  4*3  =  2  Nil  (  2-3  ) 


where  si  is  the  length  of  segment  AB,  s \  the  length  of  segment  CD,  (|)  2  and  (|)  3  are  the 
phase  changes  incurred  upon  reflection  from  the  regions  2  and  3  respectively,  ko  is  the 
free  space  propagation  constant,  and  N  is  an  integer  number.  Under  TIR  a  wave  polarized 
with  its  electric  field  parallel  to  the  interface  (TE  waves)  incurs  a  phase  shift  given  by: 


<t>  =  —  2  arctan 


N-n;kN2/kkNf-} 


1/2 


(2-4) 


For  a  wave  with  it’s  magnetic  field  polarized  parallel  to  the  interface  the  phase  shift  is 
given  by: 


4>  =  -  2  arctan 


1/2 


(2-5) 


Figure  2.2  Illustration  of  the  phase  condition  that  leads  to  the  eigenvalue  equation, 
from  reference  22 
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where  |3  is  the  propagation  constant  given  by: 

fi  =  nl  k0  COS0!  (2-6) 

Using  equations  ( 2-3  )  and  (  2-4  )  we  get  the  eigenvalue  equation  for  TE  waves: 

tanKj  =k(y +8)/(k2 -yd)  (2-7) 

where: 

K  =  (n2  kQ2  -  ft2) !  —  nxk0  sinGj  (  2-8  ) 

J  ={fi2-  n2k2Y  =  [( n2  -  n2)k2  —  K2]'  2  (  2-9  ) 

8  =  [f$2  -n2k2\  =\n2  -n2)k2  -  x2]1  (2-10) 

and  using  equations  ( 2-3  )  and  (  2-5  )  the  eigenvalue  equation  for  TM  waves  is: 

tan xd  =  n2 x {n3 y  +  n2h ) / {n2n3x 2  - n,4 y8 )  (2-11) 

Equations  (  2-7  )  and  (  2-11  )  establish  the  conditions  under  which  optical  rays  will 
propagate  through  a  slab.  These  equations  can  be  normalized  in  order  to  represent  a 
combination  of  given  parameters.  The  normalized  equations  for  TE  values  is: 


tan 


(r.'TT)  '/Ei ^)  +  V(i-feXft  +  PrTi) 

1  -  b  -  ^bJb  +  DifJv2) 


(2-12) 


and  for  TM  is: 


tan(u  VTh )  =  -b)  +  ^bjb  +  Dif2 /V2)) 

n2n2{  1  -  b )  -  n^b  [b+Dif/V2] 


(2-13) 


where  V is  the  normalized  frequency  given  by: 
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V  =  k0d  fn2  -n2  , 

b  is  the  normalized  propagation  constant  given  by: 

(lAo)2  -k22 
U  2  2  ’ 
ni  -  n2 


(2-14) 


(2-15) 


and  Dif  is  the  difference  term  given  by: 

Dif  =  k0d  -)J. ~n 22  -n32  .  (  2-16  ) 

Equations  (  2-12  )  and  (  2-13  )  will  generate  a  family  of  curves  which  represent  the 
possible  modes.  If  we  plot  these  curves  in  a  b  vs.  V  plot,  with  Dif  constant,  we  obtain 
what  is  called  a  mode  chart,  which  is  shown  on  Figure  2.3.  Using  the  mode  chart,  a 
waveguide  can  be  designed  to  guide  a  specific  number  of  modes.  We  determine  the  V 
number  required  (assuming  that  Dif  is  constant  or  determine  both  Dif  and  V)  for  the 


Figure  2.3.  Mode  chart  for  a  symmetric  slab  waveguide 
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waveguide  to  have  a  predetermined  number  of  modes.  With  this  V  (and  Dif)  number,  we 
use  equation  ( 2-14  )  to  solve  for  the  desired  parameters. 

Rectangular  Waveguides 

In  integrated  optics  more  complex  3D  structures  are  used.  A  common  structure  is  a 
rectangular  dielectric  stripe  surrounded  by  low  index  dielectric  regions.  Figure  2.4  shows 
such  a  structure.  In  here  we  see  the  dielectric  material  ni  surrounded  by  4  dielectric 
regions  along  its  sides.  A  direct  analytical  analysis  of  the  structure,  as  it  was  for  the  case 
of  the  dielectric  slab,  is  impossible.  Numerical  simulators  have  been  created  with  variable 
degrees  of  precision.  Marcatili  developed  an  approximate  analytical  solution,  for  a  wave 
guiding  structure  far  from  cut-off.  Far  from  cut-off  the  field  is  confined  mostly  in  the  core 
region  and  the  shaded  areas  can  be  neglected.  Using  Maxwell’s  equations: 


-< - c - >- 


n3 

> 

< 

X 

n5 

nx  c 

J - 

y 

>  »4 

n2 

< - b - ► 


Figure  2.4.  Rectangular  Waveguide 
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(2-17) 

(2-18) 

and  considering  time  harmonics  and  z  dependence  in  both  fields,  E  and  H ,  of  the  form: 

^l(cor-pz) 

(2-19) 

we  can  define  the  transverse 

electric  field  components  in  terms 

of  the  longitudinal 

component,  z : 

£,=-(< 

/*/)[(£  BEJBx)+  im,  dHjdy)] 

( 2-20 ) 

£,=-( 

i/Kj2)[{fi  dEjdy)—  {(Ofi0  dHjdx)] 

(2-21) 

dHjdx)- (co  n2e0  dEjdy\ 

(  2-22 ) 

»,=-(</ 

'K])[(fi  dHz/dy)+  ((On/e0  dEjdx\ 

(  2-23  ) 

where: 

Kj^jX-PT  <2-24> 

>=1,2,3, 4, 5 

and  the  subscripts  denote  the  respective  region.  The  longitudinal  components  must  satisfy 
the  wave  equation  given  by: 

{i? +  h +  K>}'v=ix’y)  = 0  (2-25) 

where  \| /2  is  the  longitudinal  electric  or  magnetic  field.  By  assuming  that  \| /2  is  separable 
and  using  equation  ( 2-19  )  we  can  write: 
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Vz  (x,y,z,t)  =  X(x)Y(y)e 


(co/-  (3z 


( 2-26 ) 


Then  the  wave  equation  can  be  written  as: 


+  +  K]  .  0 

X  dx2  Y  dy2  1 


( 2-27 ) 


by  substituting: 


Ki=Ky  +Kjy 


( 2-28 ) 


Equation  (  2-27  )  can  be  arranged  as: 


i  a2x  2  1  .  - 

IV  +  KA  +  " 


=  0 


( 2-29 ) 


A  solution  to  this  equation  is 


X(x)  =  Al  cosKjx(x+<j)x)  +  A2  sinX|x(x+(j)x  ) 

Y(y)  =  5,  cosXjy(.v+dy)  +  B2s\nK]y(y+^y) 


(  2-30 ) 


The  waveguide  is  found  to  support  a  discrete  number  of  modes  that  can  be  grouped  as 
El  and  El  modes.  E*  is  predominately  polarized  in  the  x  direction  while  Ey  is  in  the 


y  direction,  p  and  q  are  the  peaks  in  the  field  in  the  x  and  v  direction  respectively.  The 
fields  associated  with  fspxq  are  {  Ex,  Hy,  Ez,  H7,  Ey  =  0,  Hx  =  0}  and  with  Elq  are  {  Ey,  Hx, 

Ez,  Hz,  Ex  =  0,  E{y  =  0  J  .  These  modes  have  components  of  both  TE  and  TM  modes,  for 
this  reason  they  are  referred  to  as  hybrid  modes.  By  establishing  the  Ez  field  in  the  core 
for  both  modes,  is*  and  ispyq ,  in  the  form  of: 


Ez{x,y,z)  =  ^coskx(x+c|)x)cosk  (y-K|)  ) 


(2-31) 


where: 
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Kx=^lx  (2-32) 

K  y=Kly  (2-33) 

22 

and  solving  for  the  transverse  field  component  that  vanishes,  Hx,  we  can  obtain  Hz. 
Solving  equations  (  2-20  )-(  2-23  )  with  boundary  conditions  for  each  region  we  obtain 
the  eigenvalue  equations  for  E*q  as:22 


tanK,rf=  »?K,  fer. +»!?,) 

«2kx  -  y2y3 

(2-34) 

+  ^  Ky(y4+y5) 

tanK  o  =  — r - 

Ky  -  y4y5 

(  2-35  ) 

and  for  Eqq : 

,b = niKi{iy<+f^] 

n4  ns  Ky  _  n\  y4y5 

(  2-36  ) 

,a„  4  =  K  (7,  +  Y>) 

Kx2  -  y2y3 

(  2-37  ) 

22 

where: 

KX  +  Ky  =  =  Hi  ~  fi 

(2-38) 

Y2=k-n22)k02  -  k2]'/2 

(  2-39  ) 

y3  =  [(«i2-«32k2  -  kx2]1/2 

(  2-40 ) 

y4  =  [k2-«42k2  -  Ky]1/2 

(2-41) 

y5  =  k-»52k2  -<T- 

( 2-42 ) 
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The  solutions  to  equations  (  2-34  )  to  (  2-37  )  dictate  which  and  how  many  modes  will 
be  present  in  a  rectangular  waveguide.  It  is  possible  to  normalize  these  equations  to  better 
present  the  information  contained  in  them.  For  the  Ey  case  the  normalized  equations  are: 


U  =  p%  -  tan"1 


( _  \ 


u„ 


uz 

y  y 


-  tan"' 


V  1  ) 


u„ 


V  -U 

y  y 


^V2{\- B)-r2U;  =  qn  -  tan" 


If 


V2{l-B)-r2U 


V2B  +  r2  U2 


-  tan" 


V2(\-B)-r2U 


V2  B  +  r2  U2  +  Dif2 


and  for  the  E*q : 


U  =  pn  -  tan" 


U„ 


Jvf+Dip-Ul 


u 

-  tan'1 

y 

Jv2  -u2 

_  V  y  y  _ 

1p(i7sP7rc7  = 


-  tan" 


tan" 


V2 
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V2(l  -B)-r2U 


2  T  7"  2 


V  B  +  r  U 


V2(\  -  B)-  r2  U, 


V2  B  +  r2  U2  +  Dif2 


where: 


V  —  kd  n2  -  if 
Vy  =  kbjn2  -n] 
Difx  —  kd ■yj n\  -  n2 


(2-43) 


( 2-44 ) 


(2-45) 


(  2-46  ) 


(  2-47 ) 
(  2-48  ) 

(  2-49  ) 
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Dify  =  kb^n;  -  n; 

(2-50) 

II 

7\ 

(2-51) 

d  =  ^V2{\-B)-r2U2y 

(2-52) 

r  =  — ,  see  Figure  2.4 

(2-53) 

b 


B  is  the  normalized  propagation  constant,  V  is  the  normalized  frequency,  Vy  is  partial 
normalized  frequency  in  the  y  direction.  Difx  and  Dify  are  the  differential  normalized 
frequencies  in  the  x  and  y  directions.  With  these  normalized  equations  it  is  possible  to 
create  a  mode  plot  that  will  help  us  design  integrated  waveguides.  Figure  2.5  illustrates 
such  a  plot.  Although  it  looks  similar  to  the  plot  in  Figure  2.3,  it  is  only  a  2  dimensional 
slice  of  a  multidimensional  set  of  curves  governed  by  f (B,  V,  Vy,  Difx,  Dify,  r). 


B-V  Curve 


Figure  2.5.  Mode  chart  for  a  rectangular  waveguide 
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Beam  Propagation  Method 

In  the  beam  propagation  method  (BPM)  we  start  with  the  scalar  wave  equation  given 

by22: 

+  n2{x,y,z)k\f  =  0  (2-54) 

where  \|/  is  the  electric  or  magnetic  field  with  propagation  constant  k  and  traveling  in  the 
medium  with  index  of  refraction  n(x,  y,  z).  The  complexity  of  arriving  at  the  solution  is 
reduced  by  the  following  assumptions.  First,  the  variations  in  the  refractive  indices  are 
small,  and  second  the  reflected  waves  are  neglected.  The  BPM  was  conceived  on  the 
realization  that  as  optical  beam  travels  through  an  inhomogeneous  medium  it  suffers  from 
diffraction  and  from  phase  changes.  These,  although  acting  continuous,  can  be  modeled 
as  discrete  effects,  with  a  propagating  region  followed  by  a  generalized  lens.  The  lens 
would  add  the  phase  and  diffraction  effects  that  the  optical  beam  would  encounter  by 
crossing  the  particular  region.  By  assuming  that  the  wave  travels  in  the  z  direction,  the 
trajectory  can  be  divided  into  segments  of  length  A z,  each  with  a  constant  index  of 
refraction  n.  By  decomposing  the  waves  in  the  medium  between  the  segments  into  a 
spectrum  of  plane  waves  we  get:22 

\\f(x,y,Az)  =  -  *  [  U(K,a)exp[-i(ia;-i-ay-i-^Az)]dK  da  (2-55) 

(27t )  J . 

For  this  equation  the  parameters  K,  a,  and  (3  are  the  propagation  constants  in  the  x,  y,  and 
z  directions.  Since  (3  is  not  arbitrary  the  integration  does  not  cover  it.22  By  substituting 
equation  (  2-55  )  into  equation  (  2-54  )  we  get22 
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where  n(x,  y,  z)  is  replaced  by  n  .  In  most  guided  wave  problems,  the  fastest  variation  in 
the  field  \| /  is  due  to  the  phase  variation  along  the  direction  of  propagation.  It  is 
beneficial  to  factor  this  rapid  variation  out  of  the  problem.  This  is  accomplished  by 
factoring  equation  ( 2-56  )  as  follows: 


fi  —  nk  - 


(2-57) 


with  this  equation  we  can  now  separate  the  rapid  variation,  given  by  nk,  from  the  more 

2  2  22 

slow  variation  given  by  k  and  o  .We  can  now  write  \| /  as: 

\|/  (x,  y,  Az)  =  F  e'i5kAz  (2-58) 


22 

where: 


2  2 

K  +G 

exp  l - ,  Az  e 

nk  +  Jr^k2  -K 2  —<7 2 


-i(icx  +  ay) 


(2-59 


Not  only  has  the  rapid  variation  term  been  factored  out  but  the  range  of  integration  has 
also  been  changed.  For  large  values  of  K  and  a,  (3  becomes  imaginary.  This  would  mean 
that  the  ray  would  not  only  include  plane  waves  that  propagate  in  the  z-direction  but  also 
waves  traveling  at  large  angles  from  the  z-axis  and  evanescent  waves  that  decay 
exponentially  in  the  z-direction.22  The  appearance  of  these  evanescent  waves  in  the 
calculations  is  undesirable.  The  BPM  requires  that  the  propagation  process  can  be  split 
into  a  free  space  propagation  step,  in  a  homogeneous  media,  followed  by  a  lens  step,  to 
correct  the  phase  in  the  step.  The  lens  step  requires  that  the  rays  associated  with  the  wave 
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22 

travel  essentially  parallel  to  the  z-axis.  This  assumption  is  violated  when  the  ray  angles 
become  large  and  turn  imaginary  so  that  evanescent  waves  appear.  With  these  changes  we 
have  lost  some  flexibility  and  gained  in  computational  efficiency.  Factoring  of  the  rapid 
phase  variation  enables  us  to  represent  numerically  the  optical  ray  in  a  grid  much  coarser 
than  its  wavelength.  The  computational  simplification  achieved  by  the  BPM  has  its 
price.  As  mentioned  above,  the  fields  are  constrained  to  propagate  primarily  along  the  z- 
axis,  so  that  the  waves  at  large  angles  and  especially  their  evanescent  waves  do  not  make 
significant  contribution  to  the  calculations.  The  rate  of  change  of  the  refraction  index 
must  be  small,  thus  be  a  weakly  guiding  structure.  The  backward  traveling  waves  should 
not  be  necessary  in  the  behavior  of  the  device.  Fields  that  have  a  complicated 
superposition  of  phase  variation  may  not  be  accurately  modeled  if  the  phase  variation  is 
critical  to  device  behavior.  Some  of  these  limitations  have  been  relaxed  or  eliminated 
through  various  techniques,  see  reference  [23]. 

Equation  (  2-59  )  is  the  propagation  step  of  the  BPM.  For  the  lens  step  we  use  the 
equation:22 

G'  =  exp {-  i  n(x,  y,z)k  A z}  (  2-60  ) 

which  describes  a  wave  traveling  parallel  to  the  z-axis.  Since  equation  (  2-59  )  already 
includes  the  phase  shift  exp(i  n  k),  this  should  be  excluded  from  equation  (  2-60  )  giving 
us: 

G  =  exp{-  i  \n{x,  y,  z) -  n]k  A z}  (  2-61  ) 


Multiplying  equations  (  2-59  )  and  ( 2-61  )  we  have: 
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V  -FG  (2-62) 

That  is  the  mathematical  description  for  the  BPM.  Each  step  is  comprised  of  using  the 
result  from  the  previous  step  multiplied  by  FG.  Equation  (  2-59  )  can  be  shown  to  be  a 
two  nested  Fourier  transform.  These  continuous  Fourier  transform  can  be  converted  to 
discrete  Fourier  transform  to  take  advantage  of  FFT  algorithms,  see  reference  [22] 


CHAPTER  III 


ELECTRO-OPTIC  EFFECT 


When  light  transverses  an  isotropic  media  it  will  encounter  the  same  properties  in  all 
directions.  This  means  that  the  induced  polarization  will  be  parallel  to  the  electric  field 
and  will  be  related  to  it  by  a  scalar  factor.  For  an  anisotropic  media,  the  induced 
polarization  on  the  media  depends  on  the  direction  and  magnitude  of  an  applied  electric 
field.14,26  This  implies  that  the  speed  of  propagation  of  light  in  an  anisotropic  media  is  a 
strong  function  of  the  propagation  direction  and  the  polarization  of  light.14  Crystals  are 
anisotropic  in  general.  Those  that  are  isotropic  belong  to  a  class  known  as  cubic  crystals. 
Figure  3.1  shows  a  mechanical  depiction  of  an  atom  in  a  crystal  lattice.  In  this  figure  the 


Figure  3.1  Spring  model  of  Crystal  bonding,  from  reference  [24,  25] 
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springs  and  their  constant  depict  the  bonds  and  their  strength.  Using  this  mechanical 
model  it  is  easily  seen  that  the  direction  and  polarization  of  the  incident  light  on  the 
crystal  will  affect  its  propagation  speed,  i.e.  index  of  refraction. 

To  have  a  quantitative  look  at  light  traveling  in  a  crystal  we  will  start  with  Maxwell’s 
equation  for  a  general  media: 


VxE  - 


dB 

dt 


V 


^  dD  - 

dt 


V  D  =  p 

V  B  =  0 


(3-1) 

(3-2) 

(3-3) 

(3-4) 


where  J  is  the  electric  charge  density,  B  is  the  magnetic  flux  density,  D  is  the  electric 
flux  density,  and  p  is  the  charge  density.  In  optics  most  of  the  materials  used,  and  the 
ones  of  interest  in  this  research,  are  non-magnetic.  For  a  non-magnetic  material  the 
following  relationship  holds: 

B  =  [i0H  (3-5) 

Using  the  above  equations  we  can  arrive  at  the  general  wave  equation  for  non-magnetic 
materials: 


Vx 


(vxi?)  =  - 


-Ho 


Vx 


dH 

~9t~ 


(3-6) 


For  vectorial  fields,  equation  ( 2-19  )  changes  to: 

ei(cot-kr) 


(3-7) 
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With  the  vectors  E,  D,  and  H  having  the  sinusoidal  dependence  of  equation  (  3-7  )  and 
using  equation  (  3-2  )  we  get: 

Vx(vx£)=  ju0(i>2  D  (3-8) 

The  vector  k  can  also  be  expressed  as  k  =  k k' where  «  is  the  unit  vector  and  k  =  \k  |.14 
From  here  we  can  express  equation  (  3-8  )  as:14 

Vx(vx£)=-H2[i'x  (k' x  2?)]  =  - k2  [ k'(k '  -e)-e\  (  3-9  ) 

From  equations  (  3-8  )  and  (  3-9  )  we  get: 

-k-^[E-kik’E)]=D  (3-10) 

For  a  non  absorbing  media  equation  (  3-10  )  has  the  form:14 

D  =  eE  (3-11) 

where  e  is  the  dielectric  tensor  which  has  the  form:14,26 

ex  0  0~ 

e  =  0  £y  0  (3-12) 

.0  0  £z_ 

£x,  £y,  and  £z  are  the  non-negative  principal  dielectric  constants  of  £  .  For  isotropic  (cubic) 
crystals  these  dielectric  constants  are  identical.  When  only  two  of  the  constants  are 
identical  the  crystal  is  called  uniaxial,  or  bireffingent.  When  all  three  are  different  then 
the  crystal  is  called  a  biaxial  crystal.  The  majority  of  the  crystals  used  in  electro-optics  are 
of  the  uniaxial  kind.14 

The  energy  density  stored  in  an  anisotropic  crystal  is  given  by:15 
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Ue  =—D  ■  E 
e  2 


Inserting  equation  (3-11 )  and  (  3-12  )  we  get: 


D2  D, 


D: 


2U  „  =  ^  +  ^  +  ^ 


(3-13) 


(3-14) 


if  we  let 


r  = 


D 

12UX 


(3-15) 


and  by  using  the  definition  of  the  principle  indices  of  refraction,  for  a  non-absorbing 
media: 


Figure  3.2  Index  ellipsoid  of  a  uniaxial  crystal,  from  reference  [25,26] 
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i  =x,y,z 


(3-16) 


we  have: 


2  2  2 

x  y  z 

- H - H - —  1 

2  2  2 


(3-17) 


which  describes  a  general  ellipse  with  the  general  axes  parallel  to  the  x,  y,  and  z  directions 
whose  respective  length  is  2 nx,  2 ny,  and  2 nz.  In  a  uniaxial  crystal  we  can  define  the 
principle  indices  of  refraction  as: 

Hq  /lx?  Wy  ,  Hq  Yl%  (  3“  18  ) 

so  that  now  Equation  (  3-17  )  becomes: 


x2  y2  z2 

~T  +  —  +  —  -1 
n  n  n 

o  o  e 


(3-19) 


Often  in  such  a  setup,  the  z-axis  is  called  the  optical  axis.14  To  find  the  polarization 
directions  and  the  phase  velocity  of  a  given  direction  of  propagation  follow  the  following 
procedure25’26,  see  Figure  3.2: 

•  Find  the  intersection  ellipse  between  a  plane  through  the  origin,  that  is  normal 
to  the  direction  of  propagation  k ,  and  the  index  ellipsoid,  Equation  (  3-19  ). 

•  The  two  axes  of  the  intersection  ellipse  are  equal  in  length  to  2n\  and  2m 

•  These  axes  are  parallel,  respectively  to  the  directions  of  the  Dt ,  and  D2 
vectors  of  the  two  allowed  solutions 

n\  and  «2  corresponds  to  the  two  indices  of  refraction  for  the  propagation  vector  k . 
Dx  and  D2  are  the  allowed  solutions  and  the  ellipse  axes  are  parallel  to  D,  and  D2 ,25  The 
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propagation  vector  k  makes  an  angle  9  with  the  z  (optic)  axis.  Due  to  circular  symmetry 

there  will  be  no  loss  of  generality  if  we  let  x  =  0,  and  thus,  let  the  projection  of  k 
coincide  with  the  y-axis  on  the  x-y  plane.  The  allowed  polarization  directions  will  be  the 
axes  of  the  intersection  ellipse  of  the  plate  normal  to  k  with  the  index  ellipsoid,  see 
Figure  3.2.  The  refraction  indices  of  each  polarization  will  be  given  by  n0  (on  the  x-axis) 
and  ne(Q  )  on  the  y-z  plane.  Using  Figure  3.2  we  then  obtain  the  following  relations: 

y  =  -ne(0)cos(0) 

z  =  n.(e  )sin(8)  <3'20) 


Inserting  them  into  Equation  (  3-19  )  we  get: 


cos2(0)  sin2(( 


«e(e) 


(3-21) 


a  wave  polarized  along  the  n0  ( D] )  axis  is  called  the  ordinary  wave  and  that  along  the 


»e(0  )  (D2  )  is  called  the  extraordinary  wave.  Now,  with  Equation  (  3-21  )  and  given  it’s 


orientation  with  respect  to  the  crystal  axes,  we  can  determine  the  refraction  index  of  a 
wave  traveling  through  a  uniaxial  crystal.  A  uniaxial  crystal  where  ne  >  n0  is  called  a 
positive  crystal  and  when  n0  >  ne  it  is  called  a  negative  crystal.  The  index  ellipsoid 
coordinate  frame  in  general  is  expressed  as: 


z  xy  yz  zx 

- b - b - b - 

2  2  2  2 

n  n  n  n 

zz  xy  yz  zx 


( 3-22 ) 


if  we  use  the  following  substitution,  called  the  optical  impermeability  tensor: 
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" 


( 3-23  ) 


where  i,j  =  1,  2,  3  corresponds  to  the  principal  axes  x,  v,  z.  Using  the  summation  over 


repeated  indices  convention  to  rewrite  Equation  (  3-22  )  we  get: 


15,25 


Ei,  T  A  = 1 


(3-24) 


where  x\  =x,X2=  y,  etc. 

The  electro-optic  effect  refers  to  the  change  in  optical  dielectric  properties  of  the  crystal 
induced  by  an  electric  field.  According  to  quantum  theory  of  solids,  the  optical 
impermeability  tensor,  r|ij,  depends  on  the  distribution  of  solids  in  the  crystal.15  This 
distribution  can  be  changed  by  the  application  of  an  electric  field,  which  can  also  deform 
the  ion  lattice.15  The  electro-optic  effect  can  be  defined  as  a  power  series  expansion  of  T|y 

^lij(^)_ ^lij(O)  =  ^lij  =  ^jk-^k  T  Sijkl  + ...  (  3-25  ) 


where  are  the  linear  electro-optic  (Pockels)  coefficients  and  are  the  quadratic 
electro-optic  (Kerr)  coefficients.  For  this  dissertation  work  we  will  only  be  concerned 
with  the  linear  electro-optic  effect.  The  tensor  is  given  by: 


rm 

ril2 

ril3 

r221 

r222 

r223 

r331 

r332 

r333 

r231 

r232 

r233 

r321 

r322 

r323 

ri31 

ri32 

ri33 

r311 

r312 

r313 

ri21 

ri22 

ri23 

r211 

r2\2 

r213 

( 3-26 ) 
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If  the  medium  is  lossless  and  the  frequencies  of  the  electric  fields  are  much  less  than  the 
optical  fields  then  the  indices  i  and  j  can  be  permuted.15,25  This  enables  the  use  of  the 
following  contracted  notation: 


1  =(H) 

2  =  (22) 

3  =  (33) 

4  =  (23)  =  (32) 

5  =  (13)  =  (31) 

6  =  (12)  =  (21) 


( 3-27  ) 


using  the  contracted  notation  we  can  write  the  linear  electro-optic  coefficients  as: 


r,  = 


'21 


'31 


'41 


'51 


'61 


'12 


'22 


'32 


'42 


'52 


'62 


'23 


'33 


'43 


'53 


63. 


and  the  optical  indicatrix  (index  ellipsoid)  as: 

Hj  x2  +r|2  y2  +r|3  z2  +r|4  yz  +r|5  xz  +r|6  xy  - 1 


( 3-28 ) 


( 3-29 ) 


One  should  note  that  the  resultant  elements  of  the  matrix  r±  do  not  have  the  same 
transformation  or  multiplication  properties  as  a  tensor  would.  Additionally  we  have 
reduced  the  number  of  elements  in  from  27  to  18.  The  electro-optic  effect  is  now 
given  by: 

Vt1i  =  r±  (  3-30  ) 


and  the  index  ellipsoid  in  the  presence  of  an  electric  field  is  given  by: 
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which  reverts  to  Equation  (  3-17  )  by  removing  the  electric  field.  The  principle  axes  of  the 
new  ellipsoid  do  not  coincide  with  the  unperturbed  axes  (x,  y,  z).15  This  problem  can  be 
solved  by  finding  a  new  coordinate  system  in  which  the  principle  axes  of  the  ellipsoid 
coincide  with  the  axes  of  the  new  system.  This  can  be  done  by  finding  the  eigenvalues  of 
the  ri-matrix.  In  contracted  notation  it  is  written  as: 

\  lie  V 

n 6  ^2  11 4  •  (3-32) 

Jls  114  IT. 


The  electro-optic  effect  can  be  used  to  modulate  an  optical  wave  moving  through  the 
crystal.  To  do  so,  the  crystal’s  optical  indicatrix  along  with  its  electro-optical  properties, 
given  by  the  Pockel’s  coefficients,  are  used.  With  this,  one  can  establish  the  needed 
optical  wave’s  polarization  and  trajectory  through  the  crystal  axis  for  a  desired  effect, 
when  an  electric  field  is  applied.  Several  devices  implementing  the  electro-optic  effect 


were  briefly  discussed  in  Chapter  I. 


CHAPTER  IV 


DIRECTIONAL  COUPLER 


The  directional  couplers  are  versatile  devices.  They  can  be  used  as  modulators,  filters, 
routing  switches,  and  power  splitters.  In  its  most  simple  form,  the  directional  coupler 
consists  of  two  parallel  waveguides  placed  close  to  each  other.  Light  is  launched  into  one 
of  the  waveguides  and  its  optical  energy  is  transferred  synchronously  back  and  forth 
between  the  guides  by  optical  tunneling.  To  create  a  power  splitter  it  is  only  necessary  to 
limit  the  interaction  region  (the  region  where  the  two  waveguides  are  close  to  each  other) 
to  a  length  that  would  establish  the  appropriate  power  ratio.  In  order  for  this  energy 
transfer  to  actually  occur  the  phase  velocity  (i.e.  propagation  constant)  of  the  guided 
modes  between  the  guides  must  be  equal.  The  simplest  way  to  accomplish  this  is  to  have 
two  identical  waveguides.  But  two  dissimilar  waveguides  could  have  guided  modes  that 
are  “accidentally”  phased  matched.  This  occurs  only  for  a  limited  number  of  frequencies. 
Thus,  its  ability  to  act  as  a  band  pass  filter. 

The  mechanism  for  power  transfer  is  described  by  the  couple  mode  theory.  Using  this 
theory,  an  accurate  description  is  obtained  if  both  single  mode  waveguides  are  seen  as  a 
single  two-mode  structure.  A  better  understanding  is  obtained  if  the  separation  between 
the  guides  is  reduced  to  zero.  Now  the  structure  can  support  an  even  and  an  odd  mode 
with  slightly  different  propagation  constants,  see  Figure  4.1.  These  two  modes  are 
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Figure  4. 1  Even  and  odd  mode  of  a  Directional  coupler 


orthogonal  to  each  other  along  the  structure  and,  thus,  they  do  not  exchange  energy.  This 
condition  prevails  as  we  separate  the  waveguides.  The  propagation  constant  for  each 
mode  gets  closer  to  each  other,  as  we  separate  the  guides,  and  each  mode  approaches  the 
shape  of  an  isolated  waveguide,  see  Figure  4.2  a).22 

The  directional  coupler  can  now  be  explained  as  follows.  When  light  is  launched  into 
the  directional  coupler’s  interaction  region  it  will  do  so  with  a  field  distribution  that 
resembles  the  mode  of  that  waveguide,  see  Figure  4.2.  This  mode  can  be  expressed  as  the 


Figure  4.2  a)  mode’s  field  pattern  for  the  two  separated  waveguides,  b)  sum  of  the  odd 
and  even  fields 
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sum  of  the  even  and  odd  mode  of  the  compound  structure.  Since  the  propagation 
constants  for  the  even  and  odd  modes,  (3e  and  (30  respectively,  are  slightly  different;  these 
superposition  fields  will  change  their  relative  phase  as  they  travel  along  the  directional 
coupler.  The  length  required  to  transfer  completely  the  optical  energy  from  one 
waveguide  to  the  other  is  given  by: 


m  Tt 

A- A 


(4-1) 


The  relative  phase  shift  change  can  be  expressed  as  a  pair  of  differential  equations. 
Instead  of  expressing  these  equations  in  terms  of  the  parameters  of  the  even  and  odd 
modes,  it  would  be  more  appropriate  to  express  them  in  terms  of  the  parameters  of  the 
modes  of  the  individual  waveguides,  since  they  can  be  express  as  a  superposition  of  the 


even  and  odd  modes.  These  equations  are:18,22 

/  ^  =  iA  A<>(z)  ik014(z) 
a  z 

(4-2) 

=  hMi(z)  1kioA(z) 

n  v 

(4-3) 

where  A o  and  A\  are  the  complex  field  amplitudes,  and  po  and  Pi  are  the  propagation 
constants  of  the  waveguide  modes.  The  parameters  k0i  and  Kio  are  the  coupling 
coefficients  between  the  waveguide  modes.  By  convention  the  modes  have  been 
normalized  to  carry  one  unit  of  power.18  Thus  the  power  in  mode  m  is  given  by: 

Pm  (Z)  =  Am  (“  )  Am  (Z)  =  \Am  (Zf 


(4-4) 
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In  a  real  guide,  scattering  and  absorption  are  always  present.  For  this  reason  the 
propagation  constant  takes  on  a  complex  nature.  If  we  consider  two  identical  waveguides 

1  R 

with  loss  coefficient  a  so  that: 


i 

< 

ii 

(4-5) 

o 

y 

ii 

o 

y^ 

ii 

o 

y 

(4-6) 

If  we  assume  that  light  is  coupled  into  waveguide  0  at  point  z 

=  0,  then  we  have  the 

1  R 

boundary  conditions: 

zf0(0)  =  l  and  ^(0)  =  0 

(4-7) 

The  solution  to  the  waveguide  coupling  then  becomes: 

zt0(z)  =  cos(Kcz)elPz 

(4-8) 

Ax(z)  =  -isin(Kc  z)el(5z 

(4-9) 

Using  equation  (  4-4  )  we  get  the  power  flow  in  each  equation  as: 

P0=  A0{z)Al{z)  =  cos2  {Kcz)eaz 

(4-10) 

Pl=Al{z)Al(z)  =  sin2(Kcz)eaz 

(4-11) 

Equations  (  4-10  )  and  (4-11  )  show  that  the  power  is  actually  being  transferred  back  and 
forth  between  the  waveguides.  These  equations  also  show  that  there  is  a  90°  phase  lag  in 
guide  1  with  respect  to  guide  0.  This  is  produced  by  the  mechanism  that  creates  the 
coherent  transfer  of  energy  between  the  waveguides.  The  field  in  the  driving  guide 
causes  a  polarization  in  the  dielectric  material  which  is  in  phase  with  it,  and  which 
extends  in  the  region  between  the  guides  because  of  the  mode  tail.  This  polarization  then 

1  R 

acts  to  generate  energy  in  the  mode  of  the  driven  guide. 
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Given  equations  (  4-10  )  and  (4-11  ),  the  length  needed  for  complete  power  transfer  is 
also  given  by: 


L 


7t  mK 

- + - 

2k„  k„ 


(4-12) 


where  m  =  0,  1,  2,...  The  coupling  coefficient,  Kc,  is  a  strong  function  of  the  mode  tail’s 
shape  in  the  guides.18  For  well-confined  modes  in  which  the  evanescent  waves 
contribution  can  be  neglected,  the  coupling  coefficient  is  given  by  Marcatili’s  theory 


as: 


.18,27 


K  = 


2Ky2y4  e'T4C 

My4  +Kv2) 


(4-13) 


where  b  is  the  channel  width,  c  is  the  separation,  Ky  and  |3  are  the  propagation  constants  in 
the  y  and  z  directions,  respectively,  and  y4  is  the  extinction  coefficient  in  the  y  direction 
( see  Equation  (  2-41  )  ).  In  real  situations  it  may  be  difficult  to  fabricate  two  identical 

waveguides.  If  the  difference  in  their  propagation  constant  (A|3  =  |pi  -  P2I )  is  kept  small 

1  8 

then  the  power  distribution  in  the  waveguides  is  given  by 


P0  =cos2(g  z)e~az  + 


Afi 


V  ^  J 


sin2(gz) 


P{  =%sin2(gz)e" 

9 


(4-14) 

(4-15) 


where 


18 


9  =K2  + 


V  2  J 


(4-16) 
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As  equations  (4-14)-(4-16)  shows  power  will  still  be  transferred  between  the  guides. 
This  transfer  will  be  incomplete  since  equation  (  4-14  )  has  no  zero. 

Active  Directional  Coupler 

From  Chapter  III,  if  the  directional  coupler  is  made  from  an  electro-optic  crystal  then  its 
behavior  can  be  controlled  by  an  applied  electric  field.  Starting  with  a  directional  coupler 
with  two  identical  branches,  its  behavior  is  dictated  by  Equations  (  4-10  )  and  (4-11  ). 
The  condition  for  total  transfer  of  power  is  given  by  Equation  (  4-12  ).  If  a  modulating 
signal  voltage  is  applied,  as  shown  in  Figure  1.5,  a  variation  in  the  index  of  refraction  of 
the  waveguide  will  be  introduced.  This  index  variation  creates  a  mismatch  between  the 
waveguides  in  the  form  of  a  propagation  constant  difference,  A(3.  Under  the  applied 

electric  field,  Equations  (4-14)-(4-16)  govern  the  behavior  of  the  directional  coupler, 

1  & 

and  the  condition  for  total  transfer  of  power  in  this  case  is  given  by: 

gL-n  +  mn  (4-17) 

where  m  =  1,  2,  3,  ...  By  substituting  equation  (  4-12  )  into  (  4-17  )  it  can  be  shown  that 
the  A(3  required  for  100%  modulation  is:18 

{Afi)L  =  yf?>Tl  (4-18) 

The  effective  index  of  refraction  in  a  guide  is  given  by: 18 


n 


g 


(4-19) 


were  k  is  the  propagation  constant.  Thus  the  change  in  the  effective  index  needed  for 
100%  modulation  is  given  by:18 
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An 


g 


a/3  7T 
kL 


(  4-20  ) 


and  Ans  is  a  function  of  the  electric  field  applied.  From  Equation  (  4-20  )  we  can  see  that 
the  longer  the  device  the  smaller  Ang  needs  to  be.  But  there  is  a  practical  limit  in  the 
length  of  the  device.  For  integrated  devices,  the  smaller  the  size  the  greater  the  yield.  The 
increase  in  yield  is  due  to  not  only  to  the  number  of  devices  that  can  be  placed  in  a  given 
area,  but  also  to  the  reduced  probability  of  defects  in  the  device.  An  approach  to  reduce 
the  size  of  the  directional  coupler  without  reducing  A ng  is  to  have  a  push  pull 
configuration  as  illustrated  in  Figure  4.3. 18,19  In  this  type  of  configuration  the  index  of 
refraction  is  pushed  up  in  one  waveguide  and  pulled  down  in  the  other.  Other  types  of 
configurations  exist  whose  use  depends  on  the  application  requirements.18,19 

Taking  a  second  glance  at  Equations  ( 4-14  )  to  (  4-16  )  it  is  noticed  that  they  resembled 
a  sine  (  =  (sin(x)/x)  )  function.  Rearranging  then  Equations  (  4-14  )  and  (  4-15  )  in  terms 


Figure  4.3  Push  pull  directional  coupler 
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of  A(3  L  we  get: 
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P0  =  cos" 


M2  + 


f  NIL}1 

2  y 


+ 


f  N}}2 

2  y 


sin 


(K  cZ)2  + 


fm} 

2  y 


2  A 


sin 


p,=(k,l)2 


MF  + 

A 


fyuN2 
2  y 


(kl)2  + 


2  y 


M)2  + 


A/iZ 

v  2  y 


(4-21) 


( 4-22 ) 


Setting  Kc L  =  %/2,  see  Equation  (  4-10  ),  and  plotting  the  equation  we  get  the  typical 
switching  characteristics  of  an  active  directional  coupler.19  As  it  can  be  seen  in  Figure 
4.4,  when  A|3Z  =  0  all  the  power  comes  out  of  guide  1.  As  A|3L  increases  (  and  thus  the 
electric  field  )  the  directional  coupler  switches  to  the  off  state  (A(3Z  =7tV 3 ).  If  A(3L  is 


Figure  4.4  Switching  a  directional  coupler 
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increased  further  the  coupling  between  the  guides  will  be  destroyed.  From  Figure  4.4  we 
can  deduce  that  the  off  state  is  electrically  adjustable.  If  in  the  initial  state  Kc L  is  not  made 
nil,  by  a  fabrication  error,  then  the  two  switching  states  could  not  be  obtained  form  the 
configuration  on  Figure  1.5  or  Figure  4.3. 19  To  solve  this  problem  there  exist  a  different 
configuration,  known  as  A(3  reversal,  in  which  an  extra  set  of  electrodes  are  added  to  the 
structure,  but  with  inverse  polarity,  see  references  18  and  19.  Each  set  of  electrodes  is 
independently  adjusted  and  with  LI 2  length. 

Symmetrical  Directional  Coupler 

At  the  beginning  of  the  chapter  a  zero-gap  directional  coupler  was  envisioned  to  explain 
the  behavior  of  these  devices  in  general.  Far  from  being  a  theoretical  curiosity  these  are 
viable  devices.  In  the  literature  at  least  three  names,  symmetric  directional  couplers  , 
BOA  switches  ,  and  zero-gap  directional  coupler  have  referred  to  these  devices.  The 
power  flow  equations  that  govern  the  symmetric  directional  couplers  are  the  same  as 
those  for  the  identical  guides  directional  coupler,  Equations  (  4-10  )  and  (4-11  ).  The 
value  of  Kc,  as  it  was  for  the  former  case,  is  given  by: 

kc=|(A-A)  (4-23) 

also  it  is  assumed  that  the  loss  coefficient  for  each  waveguide  is  identical.  If  a  change  5n 
is  introduced,  by  way  of  the  electro-optic  effect,  in  the  waveguide  index  the  mode  shapes 
are  altered  in  a  symmetric  way,  they  will  remain  orthogonal.  With  this  change  the 
perturbed  propagation  constants  will  differ  from  the  unperturbed  ones,  generally  by 
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different  amounts.  Therefore  equations  (  4-10  )  and  (4-11  )  would  still  hold  for  the 
perturbed  directional  coupler  but  with  a  different  coupling  coefficient  which  is  given  by: 

<=Kc+Kceo  (4-24) 

where  Kceo  is  the  change  in  the  coupling  coefficient  due  to  a  perturbation,  the  electro-optic 

effect  in  this  case.  Starting  with  a  directional  coupler  with  complete  power  transfer,  the 
necessary  change  to  switch  symmetrically  a  directional  coupler  is  given  by: 


(4-25) 


(4-26) 


( 4-27  ) 


The  term  d(fie  - fi0)/dn  depends  strongly  on  the  separation  of  the  waveguides  and 

2 1 

increases  as  the  separation  decreases.  This  term  has  been  shown  to  increase  by  an  order 
of  magnitude  from  an  ordinary  directional  coupler.  This  makes  it  possible  to  switch  a 
zero  gap  coupler  with  a  symmetric  change  of  5n  with  reasonable  voltages.  In  the  case  of  a 

directional  coupler  with  a  waveguide  separation,  a  symmetrical  change  of  8n  would 

2 1 

require  excessively  high  voltages. 
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Since  the  only  thing  that  changes  in  the  power  flow  equations  is  the  coupling 
coefficient,  it  is  now  possible  for  the  directional  coupler  to  achieve  both  states,  on  or  off, 
electrically,  and  with  only  one  controlling  voltage.  This  makes  the  devices  able  to 
compensate  for  fabrication  errors  by  applying  a  bias  voltage.  This  would  be  a  much 
simpler  arrangement  than  using  the  A(3  reversal  configuration  of  a  traditional  coupler. 

Crosstalk 

Crosstalk  in  a  directional  coupler  is  defined  as  the  ratio  of  light  power  in  the  unwanted 

AA 

output  port  to  the  power  in  the  desired  output  port.  Another  measure  of  crosstalk  is  the 
extinction  ratio.  The  extinction  ratio  is  defined  as  the  ratio  of  the  optical  power  in  the 
“on”  state  to  the  optical  power  in  the  “off’  state,  in  the  same  guide  branch.  In  a 
directional  coupler  it  is  desired  that  all  the  power  be  coupled  to  the  other  waveguide  in 
the  “off’  state  and  none  in  the  on  state.  But  in  reality  this  is  not  the  case,  there  is  always  a 
small  amount  of  unwanted  residual  power  in  the  off  waveguide.  The  source  of  the 

AO 

crosstalk  is  due  to  different  mechanisms.  Among  them  are  unequal  mode  excitations  , 
non-optimal  coupling  length19,  absorption  loss30,  waveguides  asymmetry28,  parasitic 
coupling31,  etc. 

Physical  asymmetries  in  the  structure,  which  arise  during  fabrication,  will  lead  to 

AO 

crosstalk  due  to  coupling  between  the  modes.  These  asymmetries  manifest  themselves 
as  differences  in  the  width,  index  of  refraction,  and  curve  shape.  In  a  perfectly  symmetric 
directional  coupler,  the  even  mode  field  distribution  \|/e  at  any  cross  section  (  z  =  z\,  see 
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Figure  4.1)  is  orthogonal  to  the  odd  mode  field  distribution  \|/0  at  any  other  cross  section 
( z  =  Z2  =  z\  +  5z).28,29  In  other  words:28’29 

}¥e(f.ziVo(f.z2)  =  0  for  any  zvz2  (4-29) 

In  the  presence  of  asymmetries  there  is  no  longer  an  even  and  an  odd  mode  for  the  total 
structure.29  When  the  gap  between  the  waveguides  is  zero  or  small,  a  slight  structure 
asymmetry  will  produce  only  a  small  modification  of  the  electric  field  profiles.28,29  These 
deviations  increase  as  the  separation  between  the  waveguides  is  increased.  The  mode  in 
each  waveguide  will  start  to  resemble  the  mode  of  an  isolated  waveguide,  faster  than  with 
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Figure  4.5  Effects  of  asymmetries  in  a  directional  coupler.  Eigenmode  profiles  for 
various  fabrication  differences.  From  reference  28 
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98  9Q 

no  asymmetries  present.  ’  Thus  the  smaller  the  guide  separation  the  smaller  the 
influence  of  any  given  asymmetry.  ’  When  there  is  a  separation  between  the 
waveguides,  it  has  been  shown  that  even  small  variations  in  the  asymmetries  cause 
significant  deviations  in  the  combined  structure’s  modes.  ’  Figure  4.5  (  from 
reference  28)  shows  different  electric  field  profiles  arising  from  fabrication  asymmetries, 
index  of  refraction  and  width  differences.  Additionally,  the  figure  shows  the  effect  of  the 
gap  separation  in  the  electric  field  profiles. 

Any  waveguide  discontinuity  or  irregularity  that  differentially  affects  the  mode 

98 

propagation  losses  will  have  adverse  effects  on  the  directional  coupler  performance. 
This  would  create  an  imbalance  in  the  mode’s  amplitude  that  will  manifest  as  crosstalk  in 
the  directional  coupler.  Crosstalk  for  a  lossy  material  has  been  studied  theoretically  (  see 
reference  30)  but  on  a  two,  finite  separation,  waveguide  directional  coupler.  In  that  study 
adjusting  the  coupling  lengths  could  not  compensate  the  effects  of  differential  modal 
absorption.  The  data  presented  tends  to  indicate  that  a  reduction  in  the  waveguide 
separation,  the  crosstalk  due  to  differential  absorption  would  reduce  its  contribution  to  the 
overall  crosstalk.  A  point  to  consider  is  that  the  reduction  in  the  guide  separation  also 
reduced  the  coupling  length  necessary  for  complete  power  transfer. 

Bending  Sections 

The  theory  presented  thus  far,  and  in  most  publications  and  textbooks,  deals  with  the 
parallel  coupling  section  and  neglects  the  bending  sections,  Figure  4.6.  The  bending 
sections  of  a  directional  coupler  should  be  a  slowly  converging  or  diverging  curve  to 
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minimize  radiation  losses.28,31  The  slow  convergence  of  the  bending  sections  result  in  a 
longer  effective  coupling  length  for  the  device.  From  Equation  (  4-13  ),  as  the  two  guides 
approach  each  other  the  coupling  coefficient,  Kc,  increases.  When  Kc  is  large  enough  the 
guides  will  start  accumulating  an  appreciable  amount  of  propagation  phase.28,31  A  slow 
curved  convergence  will  accumulate  much  more  propagation  phase  than  straight 
convergence,  as  demonstrated  in  reference  28.  Therefore  the  input  bending  sections  needs 
to  be  taken  into  consideration  when  designing  a  directional  coupler.  The  generalized 
propagation  phase  angle  is  given  by:31 

^(z)  =  |oKc(z)az  (4-30) 


where  in  the  Marcatili  approach  Kc  is  given  by,  see  Equation  (  4-13  ): 

2  k  2v  e"Y4c(z) 

„  (  \  ZIS  14  C 

‘  Z  “M42  +  k/) 


(4-31) 
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Another  factor  that  needs  to  be  taken  into  consideration  is  the  bending  section  angle 
0  /2,  see  Figure  4.6.  This  angle  establishes  an  optical  wave  front  tilt  that  is  often 
neglected  in  the  analysis  of  the  directional  coupler.  The  tilted  wave  fronts  result  in 
reduced  waveguide  coupling.  The  tilt  introduces  phase  differences  in  the  regions  where 
the  coupling  takes  place,  which  decreases  the  waveguide  interaction.  An  estimate  of  the 
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permissible  tilt  can  be  obtain  by,  see  Figure  4.6: 

.  (e^  x 

sm  -  <— —  (4-32) 

^2  J  4 Wn  v 

where  W  is  the  waveguide  width,  n  is  the  index  of  refraction  of  the  core  channel.  This 
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equation  is  based  on  the  assumption  of  well-confined  modes. 

By  having  curved  bending  sections  the  angle  of  the  instantaneous  slopes  of  the  curves  is 
less  than  if  the  same  section  was  replaced  with  straight  ones,  over  the  section  where 
significant  mode  interaction  takes  place  (~  2  waveguide  width  for  a  well  confined 
modes).  This  permits  the  use  of  steeper  equivalent  tilt  angles  before  the  tilt  effects 
diminish  the  waveguide  coupling.  Additionally,  the  gradual  approach  of  the  guides 
minimizes  the  differential  mode  propagation  loss. 


CHAPTER  V 


AlGaN 


The  AlxGai.xN  alloy  system  can  have  a  wurtzite  (hexagonal)  or  a  zincblende  (cubic) 
lattice  structure.  Wurtzite  is  the  preferred  structure  due  to  the  better  crystalline  quality 
that  can  be  obtained  with  it,  at  the  present  time.  ’  The  hexagonal  unit  cell,  illustrated  in 
Figure  5.1,  is  defined  by  two  coplanar  vectors,  ai  and  a2,  at  120°  to  one  another  and  a 
third  axis  c  at  a  right  angle,  see  Figure  5.1. 34  The  complete  lattice  is  built  by  shifting  the 
unit  cell  in  the  plane  defined  by  the  ai  and  a2  vectors.  The  unit  cell,  as  seen  on  Figure  5.1, 
has  a  six-fold  rotational  symmetry  along  c,  which  generates  the  hexagonal  symmetry  of 


Figure  5.1  Hexagonal  Unit  Cell  (heavy  lines) 
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the  lattice.34  The  third  axis,  of  the  hexagonal  symmetry,  which  lies  in  the  base  plane  of 
the  hexagonal  prism,  is  related  to  a\  and  (12  in  such  a  way  that  it  is  often  used  with  the 
other  two  axes.  Thus  the  Miller-Bravais  indices  of  a  plane  in  a  hexagonal  system  refer  to 
four  axes  written  as  (hkil).  Where  h,  k,  i,  and  /  are  the  reciprocal  of  the  fractional  intercept 
on  the  <*1,  «2,  «3,  and  c  axes  respectively.34 

The  wurtzite  AlxGai_xN  alloy  system  covers  a  direct  and  wide  band  gap  range.  This 
range  expands  from  the  visible,  3.4  eV,  to  the  UV,  6.2  eV,  for  x  =  0  to  1  respectively.  For 
this  reason,  the  short  wavelength  opto-electronic  applications  have  fueled  the  interest  in 
this  material  system.  Another  reason  is  its  advantages  for  high  power  electronics, 
including  the  large  bandgap  for  high  voltage  and  temperature  operations.  The  material 
parameters  of  GaN  and  AIN  are  listed  in  Table  1  along  with  other  semiconductors  for 
comparison. 

Up  to  this  date,  a  native  substrate  for  the  AlxGai_xN  material  system  has  been  hard  to 
achieve,  mainly  due  to  the  large  equilibrium  vapor  pressure  of  nitrogen  over  the  liquid 
nitride  phases.  Bulk  single  crystal  growth  has  been  limited  to  only  ~100  mm  with  good 


Table  1.  Material  Parameters  for  300  K 


Property 

GaN 

AIN 

SiC 

3C/6H 

GaAs 

Energy  gap  (eV)(direct) 

3.4L35J 

6.2L35J 

2. 2/2. 9  [10J 

1.423[10] 

Density  (g/cm3) 

6.1[35] 

3.26[35] 

5.3 1 7[36] 

Dielectric  constant  (e) 

00 

OO 

9.7[10] 

12.85[36] 

Electron  mobility 
(cm2/V-s) 

1000  bulk[35] 

135[37] 

1000/600[1°] 

8500[36] 

Hole  Mobility  (cm2/V-s) 

30[35] 

14[35] 

40[10] 

400[36] 

Saturation  Velocity  (m/s) 

2.5xl07  [37] 

1.4xl07  [37] 

2xl07  [10] 

2xl07  [10] 

Breakdown  field  (V/cm) 

>5x1 06  [37] 

4xl06  [10] 

0.4xl06  [10] 
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quality,  but  in  a  limited  quantity.39,40’41,42  For  this  reason,  alternate  substrates  have  been 
used. 

Among  the  alternative  substrates  are  sapphire,  SiC,  GaAs,  and  Si.  Sapphire  has  been  the 
most  widely  used  substrate  to  date,  mainly  due  to  its  wide  availability.43  Although  from 
the  lattice  mismatch  (~14%),44  sapphire  should  have  been  a  poor  candidate.  SiC  is  a  much 
better  candidate  for  a  substrate  because  its  lattice  mismatch  is  only  3.4%.  SiC  has  failed 
to  produce  better  performance  than  films  grown  on  sapphire,  although  it  is  improving.45 
This  has  made  sapphire  the  predominant  substrate  for  the  AlGaN  material  system.  A 
characteristic  that  sets  this  material  system  apart  from  other  semiconductors  is  the  fact 
that  in  spite  of  the  high  dislocation  density  (108-1010  cm2)44,46,  which  is  orders  of 
magnitude  higher  than  what  it  is  observed  in  other  opto-electronic  materials,  it  has 


Figure  5.2  Threading  dislocations  in  GaN  Grown  by  MOCVD  from 
reference  47 
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obtained  high  optoelectronic  performance.44  Figure  5.2  shows  a  Transmission  Electron 
Microscopy  (TEM)  view  of  a  GaN  sample  from  reference  47. 

Crystal  Growth 

There  have  been  several  approaches  in  growing  AlGaN  epitaxy  on  a  number  of 
substrates.  Among  them,  Hydride  Vapor  Phase  Epitaxy  (HVPE),  Metalorganic  Chemical 
Vapor  Deposition  (MOCVD),  and  Molecular  Beam  Epitaxy  (MBE).  Yet  it  wasn’t  until 
Yoshida  recognized  the  importance  of  a  low- temperature  buffer  layer  for  the 
improvement  of  surface  morphology  and  crystalline  quality  of  GaN,  that  devices  in  this 
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material  system  begun  to  take  off.  ’  Subsequent  studies  improved  on  this 
approach.49,50,51  The  approach  has  a  few  hundred  angstroms  of  AIN  deposited  at  a  low- 
temperature  (~600°C).  The  thickness  and  temperature  is  optimized  for  a  particular  reactor 
and  substrate  type  49,52  Subsequent  layers  are  grown  on  top  of  this  buffer  layer.  The  buffer 
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layer  relieves  the  strain  between  the  substrate  and  subsequent  AlGaN  layers.  ’  The  low 
growth  temperature  of  the  buffer  layer  causes  it  to  be  amorphous.10  This  permits  the 
buffer  to  cover  the  substrate  uniformly.10  When  the  heat  is  increased  to  normal  growth 
temperature  the  buffer  layer  crystallizes  and  provides  a  good  template  for  further  epitaxy. 
For  sapphire  substrates  the  improvements  obtained  with  the  AIN  buffer  layer  were:10 

1 .  Two  orders  of  magnitude  decrease  in  background  carrier  concentration. 

2.  A  factor  of  10  in  increased  mobility. 

3.  Band-Gap  photoluminescence  (PL)  were  two  orders  of  magnitude  more  intense. 


4.  Suppression  of  midgap  emission 
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5.  The  x-ray  diffraction  peak  was  four  times  narrower. 

Low  temperature  GaN  buffer  layers  have  also  been  used  with  no  clear  advantage  of  one 
over  the  other,  although  laser  diodes  have  been  obtained  with  GaN  buffer  layers  with 
MOCVD. 10,43 

Of  the  different  growth  techniques  the  ones  that  have  grown  the  highest  quality 
materials  have  been  MBE  and  MOCVD.  For  commercial  applications,  MOCVD  is  the 
leading  candidate  with  blue  LEDs  and  violet  LD  already  being  marketed.  ’  ’  In 
MOCVD,  trimethylgallium  (TMGa,  Ga(CH3)3)  and/or  trimethylaluminum  (TMA1, 
A1(CH3)3)  react  with  NH3  at  a  substrate  which  has  being  heated  to  ~1000°C  to  obtain  the 
basic  reaction  given  by10,38 

Ga(CH3)3(v)  +  NH3(v)  =  GaN(s)  +  3CH4(v)  (  5-l  ) 

The  details  of  this  reaction  are  not  yet  well  understood  and  the  intermediate  steps  are 


N  2  +H  2 


Figure  5.3  Dual  Flow  MOCVD,  from  reference  66 
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thought  to  be  complex.  In  a  quest  for  superior  material  growth  Nakamura  designed  a 
special  atmospheric  pressure  MOCVD  reactor  optimized  for  GaN  growth,  see  Figure 
5. 3. 66  In  this  approach  the  reactants  are  diluted  in  H2.  The  sources  are  supplied  from  the 
horizontal  inlet  and  a  vertical  subflow  then  drives  the  reactants  to  the  growing  film 

OO  /T/T 

surface.  ’  In  MOCVD,  the  high  growth  temperature  required  to  break  the  N-H  bond  of 
ammonia  precursors  creates  a  problem.  Due  to  the  thermal  mismatch  with  all  the 
available  substrates,  post  growth  cooling  introduces  significant  amounts  of  strain  and 
defects  into  the  nitride  film.10  Additional  problems  of  the  MOCVD  reactors  are  the  film 
nonuniformity  and  the  tendency  of  ammonia  to  prereact  with  the  group  III  metalorganic 
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compounds  to  form  nonvolatile  compounds.  These  problems  need  to  be  overcome  in  the 

O  O 

design  and  scale  up  of  AlGaN  deposition  system  by  MOCVD  reactors. 

Another  technique  for  growing  AlGaN  is  MBE.  MBE  is  in  essence  a  specialized  form 


HEATER  SUPPLY 


Figure  5.4  Basic  MBE  design  from  reference  54  and  55 
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of  vacuum  deposition  with  several  important  differences.55  The  growth  of  the 
semiconductor  fdm  takes  place  with  the  reaction  of  molecular  beams  of  the  constituent 
elements  with  a  crystalline  substrate  surface  held  at  a  suitable  substrate  temperature  under 
UHY  conditions.55  The  molecular  beams  are  generated  from  sources  contained  in  effusion 
cells,  which  are  aimed  at  the  heated  substrate.55  A  basic  MBE  system  is  shown  on  Figure 
5.4. 54,55  For  the  AlGaN  growth,  gallium  and  aluminum  are  supplied  by  evaporation  from 
solid  source  effusion  cells,  and  the  atomic  nitrogen  species  are  obtained  by  either  an  N2 
plasma  source  or  by  thermally  activated  decomposition  of  NH3  on  the  growing  surface.43 
An  advantage  for  MBE  grown  material  is  the  low  growth  temperature  needed.  ’  ’  This 
should  result  in  lower  thermal  stresses  upon  cooling,  less  diffusion,  and  reduced  alloy 
segregation.  Another  advantage  is  that  as  grown  p- type  GaN  films  don’t  require  any 
passivation  steps,  due  to  the  absence  of  atomic  hydrogen  species  during  growth.  Although 
most  of  the  high  quality  material  has  been  grown  by  MOCVD,  MBE  has  made  a  lot  of 
progress  toward  AlGaN  growth,  as  is  evident  by  the  realization  of  LEDs  and 
electronics.57,58,59,60 

Recently,  the  Lateral  Epitaxy  Overgrowth  (LEO)  technique  has  attracted  attention  as  a 
way  to  reduce  threading  dislocations  in  selected  areas  of  AlGaN  layers.  ’  In  this 
technique  a  mask  is  deposited  on  top  of  the  standard  dislocated  layers  and  patterned.63,64 
A  regrowth  is  done  on  the  pattern.  The  material,  which  grows  vertically  through  the 
opening  of  the  mask,  retains  the  same  dislocation  density.  But,  the  material  that  grows 
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laterally,  on  top  of  the  mask,  has  a  much  lower  dislocation  density.  ’  ’  Figure  5.5 


Figure  5.5  TEM  showing  dislocation  density  reduction  via 
LEO,  from  reference  61 

shows  the  dislocation  density  reduction  on  top  of  a  Si02  mask.  Devices  have  been  grown 
on  top  of  the  LEO  material  with  significantly  improved  performance  compared  to  devices 
grown  on  normal,  dislocated  materials.63’ 46 

Another  big  catalyst  in  the  emergence  of  nitride  devices  was  the  achievement  of  low- 
resistivity  p- type  doping  in  the  AlGaN  material  structure.  This  achievement  was 
accomplished  by  Amano  et  al.  in  1989  using  Low-Electron  Beam  Irradiation  (LEEBI) 
treatment  with  Mg  doping  on  MOVPE  grown  material.65  Though  at  that  time,  the 
acceptor  activation  ratio  was  low  requiring  a  high  Mg  incorporation  for  high  />doping 
levels.  Later  Nakamura  et  al.  were  able  to  improve  the  p- type  doping  levels  (~  10  /cm  ) 
on  MOCYD  grown  material.66  A  deficiency  of  the  LEEBI  treatment  is  that  it  only 
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activated  a  thin  surface  layer.65,66  This  was  later  overcome  by  Nakamura  et  al.  by 
discovering  that  p- type  levels  could  also  be  activated  by  a  thermal  anneal  at  700°C  in  an 
N2  atmosphere,  this  time  the  entire  Mg  doped  layer  became  p- type  with  low  resistivity. 
This  activation  is  reversible  when  the  material  is  annealed  in  NH3  atmosphere.9  This 
reversibility  pointed  to  hydrogen  as  a  main  factor  in  preventing  p- type  doping.  As  it  turns 
out,  during  growth  hydrogen  is  interstitially  incorporated  creating  an  H-Mg  acceptor 

O  O 

complex.  This  complex  causes  the  acceptor  to  be  passivated.  This  bond  has  to  be  broken 
in  order  to  achieve  p- type  material.  This  is  achieved  with  high  temperature  annealing. 
With  this  knowledge,  p- type  material  was  also  achieved  in  MBE  grown  material  by 
minimizing  the  effect  of  hydrogen.9,68 

Optical  Properties 

As  mentioned  before,  AlxGai_xN  has  a  large  bandgap  starting  at  3.4  eV.  For 


Transmission  vs.  Wavelength 


Figure  5.6  AlGaN  transmissivity 
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semiconductors  this  makes  the  material  transparent  to  wavelengths  whose  energy  is 
below  that  of  the  bandgap.  The  transmisivity  of  AlxGai_xN  for  several  values  of  x  is  given 
in  Figure  5.6.  Along  with  this,  the  bandgap  can  be  changed  by  varying  x.  The  variation  in 
x  will  also  change  the  refraction  index  of  the  material.  Control  of  the  index  of  refraction 
enables  one  to  confine  the  optical  wave,  as  mentioned  in  chapter  II.  Precise  knowledge  of 
the  index  variation  with  respect  to  x  and  the  wavelength  is  needed.  There  have  been 
several  studies  for  the  index  of  refraction  of  AlGaN.  ’  ’  ’  These  have  inconsistencies 
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among  them  regarding  the  values  of  the  refraction  index.  This  is  due  to  the  fact  that 
these  semiconductors  are  still  in  their  infancy.  The  AlGaN  wurtzite  crystal  is  uniaxial, 
therefore  it  has  two  indices  of  refraction  related  to  it,  see  chapter  III.  Bergmann  et  al.  have 
done  a  recent  study  on  both  ordinary  and  extraordinary  refractive  indices  on  AlGaN.69 
The  measurements  were  taken  for  the  wavelengths  457<  X  <  980  nm  and  the 
compositions  x  =  0.00,  0.04,  0.08,  0.11,  and  0.20  with  an  accuracy  of  ±10%.69  The 
refraction  indices  were  obtained  using  the  prism-coupled  waveguide  technique,  and  had 
an  accuracy  of  ^0.01. 69  The  data  obtained  was  fit  to  the  1st  order  Sellmeier  dispersion 
formula  given  by: 

.  ,  A.  X2 

n(X)-l  +  X2  -A2  (5‘2) 


where: 


A0(x)  =  B0  +  Bxx  +  B2  x2 
X0  (x)  =  C0  +  Cj  x 


(5-3) 

(5-4) 
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Table  2.  Coefficient  for  the  Sellmeier  dispersion  formula,  from  reference  69 


Coefficient 

n0 

ne 

Bo 

4.1416  ±0.008 

4.3076  ±  0.008 

Bi 

-1.0  ±0.2 

-1.2  ±0.2 

b2 

-4.4  ±  0.9 

-3.7  ±0.8 

Co 

187.4  ±  0.4  nm 

190.3  ±  0.4  nm 

Ci 

-121  ±  3  nm 

-121  ±4  nm 

Table  2  lists  the  coefficients  for  the  Sellmeier  formula.  The  major  source  of  problems  in 


establishing  the  accuracy  of  the  refractive  index  equations  is  not  only  the  index  data  itself 
but  also  the  accuracy  of  the  A1  mole  fraction  in  the  AlxGai_xN  material.  According  to 
Laws  et  al.,  the  determination  of  the  A1  composition  on  AlxGai_xN  is  very  difficult  and  the 
results  depend  on  the  homogeneity  of  the  samples  as  well  as  the  measurement  technique 
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employed. 

GaN  ‘s  hexagonal  crystal  symmetry  corresponds  to  the  6mm  group.  This  enables  it  to 
have  linear  electro-optic  (Pockels)  effect.  Its  contracted  electro-optic  tensor  (see  Chapter 
III)  is  given  by: 


substituting  this  equation  into 
here  for  convenience: 


f  1  r  1 

2  rnEz 

r  i  r  i 

2  rn^z 

x2+ 

l"-  J 

{n0  J 

^  0  0  r13  ^ 

0  0  rn 

0  0  r33 

0  r5i  0 

r51  0  0 

0  0  0 

V  J 

equation  (  3-3 1  )  we  get  equation  ( 


(5-5) 


1-5  ),  which  is  shown 


/+ 


2  r33^Z 


\ 

z1  Ylyz  rslEy+  2xz  r5lEx=  1 

) 


(5-6) 
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This  equation  gives  one  the  change  in  the  index  of  refraction  as  a  function  of  the  applied 
electric  field.  Long  et  al.  has  measured  and  reported  on  some  of  the  nonlinear 
coeffitients.16  These  are  rn  =  ~  0.57  ±  0.1 1  pm/V  and  r33  =  ~  1.91  ±  0.35  pm/V.  The  third 

electro-optic  coefficient  was  measured  by  Khan  et  al.  and  is  given  by  r5l  =  3.1  ±0.2 
pm/V.17  All  these  were  measured  for  the  wavelength  632.8  nm.  Not  only  does  it  have 
linear  electro-optic  but  also  second  order  non-linearities,  which  are  beyond  the  scope  of 
this  work.73’74 


CHAPTER  VI 


EXPERIMENTATION  AND  RESULTS 


The  proposed  directional  coupler  design  for  this  work  is  based  on  a  rectangular  dielectric 
waveguide  with  zero  gap  between  the  branches,  see  Chapter  IV,  and  is  based  on  the 
AlxGaj.xN  material  system.  The  waveguides  will  have  a  clad-core-clad  configuration,  as 
shown  in  Figure  6.1  a).  The  advantages  of  this  structure  are  that  the  top  and  bottom 
cladding  can  be  doped  p  and  n  while  the  core  is  insulating.  This  will  create  a  PIN 
structure  that  in  reverse  bias  will  concentrate  the  electric  field  E  in  the  core  region.  This 
would  result  in  greater  index  variation  for  a  given  voltage.  In-house  material  was  grown 
using  an  MBE  machine  on  top  of  sapphire  substrates.  For  the  in-house  MBE,  p- type 
doping  was  not  feasible.  Therefore,  instead  of  a  p-i-n  structure  an  n-i-n  structure  is 


n/p 

I 

p/n 

clad 

air 

core 

clad 

(a)  (b) 

Figure  6. 1  Proposed  structure  a)  single  waveguide  b)  interaction  region 
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chosen. 


Design 

For  the  design  of  the  optical  waveguides,  as  stated  in  the  chapters  above,  an  index  of 
refraction  variation  is  needed  for  the  confinement  of  light.  At  the  beginning  of  this  work 
the  studies  made  by  Brunner  et  al.  were  used  as  a  basis  for  the  index  of  refraction 
calculations.71  The  wavelength  selected  for  the  operation  of  the  directional  coupler  was 
632.8  nm.  This  was  based  to  the  fact  that  the  available  information  on  GaN’s  linear 
electro-optic  coefficient  was  done  at  this  wavelength  and  the  accessibility  of  HeNe 
lasers.16  The  equations  needed  to  obtain  the  index  of  refraction,  according  to  Brunner  et 
al,  are:71 


nl(hv,x)  =  C W+  A(x)yi 2-  (l  +  yT  ~  ft  -yT ) 

He  ~  Ho  +  0-02 


(6-1) 


where: 


hv 


C(x)  =  const(x) 


(6-2) 


Table  3.  Refraction  index  parameters  for  Equation  ( 

6-1) 

x  (Al  composition) 

A(x) 

C(x) 

Bandgap  (eV) 

0.00 

9.76472 

2.74493 

3.42 

0.11 

10.04172 

2.59099 

3.64 

0.20 

11.65128 

2.04179 

3.89 

0.38 

11.87005 

1.96916 

4.16 

0.50 

12.72810 

1.25847 

4.63 

0.65 

12.56483 

1.21344 

4.73 

0.79 

12.64459 

0.78265 

5.21 

0.86 

12.64948 

0.64295 

5.50 

1.00 

13.50000 

0.75000 

6.13 
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E„(x)  is  the  effective  bandgap.  A(x)  is  a  function  of  the  oscillator  strength  of  the  optical 
transition,  and  E„(x).  A(x)  and  C(x)  experimental  values  are  given  in  Table  3.  These 
numbers  were  supplied  by  Dr.  O.  Ambacher  and  are  part  of  reference  71.  With  Equation  ( 
6-1  )  we  can  calculate  the  index  variation  as  a  function  of  A1  composition  for  each  of  the 
A1  concentrations  in  Table  3.  In  Figure  6.2  the  diamonds  are  the  calculated  refraction 
indices  and  the  solid  line  is  the  trend  line.  The  trend  line  equation  is  given  by: 

na  =  3.9968  x5  -  7.0603  x4  +  3.8812  x3  -  0.9418  x2  -  0.1613  x  +  2.3024  (  6-3  ) 

With  this  index  variation  in  terms  of  A1  concentration  we  can  now  design  waveguides  in 
a  clad-core-clad  configuration. 

As  it  can  be  seen  from  Figure  6.2  the  refractive  index  decreases  as  the  A1  concentration 


0  0.2  0.4  0.6  0.8  1 


A1  concentration  (x) 


Figure  6.2  Refraction  index  variation  for  E  =  632.8  nm  from  ref  71 


66 


increases.  For  this  reason  GaN  is  chosen  for  the  core  material.  Due  to  growth  constraints 
10%  A1  concentration  was  chosen  in  the  beginning.  The  refraction  indices  given  by 
Equation  (  6-3  )  for  the  core  and  cladding  are  then 

n0  ne 

Core,  GaN  2.342  2.362 

Cladding,  Al0.iGa0.9N  2.320  2.340 

Cladding,  Al0.05Ga0.95N  2.332  2.352 

In  a  non-native  substrate  like  sapphire,  the  Al^Ga^N  material  system  grows  in  the  c-axis 
direction,  see  Figure  5.1.  This  means  that  the  electric  field  will  be  along  the  z  direction  in 
the  crystal  axes  and  simplifies  Equation  (  5-6  )  to: 


—  +  ruEz  x2  +  —  +  rnEz  y2  +  —  +  rnEz 


(6-4) 


Since  of  these  two  r33  is  the  highest  coefficient,  the  polarization  direction  was  chosen  to 
be  along  the  z-axis.  That  will  give  us 


Figure  6.3  Cross  sectional  view  of  a  directional  coupler  in  the  interacting  region 
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2  +  r33^Z  —  1  (  6-5  ) 

"e 

Which  corresponds  to  the  e-ray  polarization. 

In  order  for  the  directional  coupler  to  work,  each  of  the  waveguides  needs  to  be  single 
mode  and  the  region  were  the  two  waveguides  are  joined  needs  to  be  dual  mode,  see 
Chapter  IV.  We  can  first  start  with  the  Marcatili  model,  since,  it  can  give  us  a  rough 
estimate  of  the  waveguides  dimensions,  using  Figure  6.3  and  Equations  (  2-45  )  -  ( 2-53  ), 
with  little  computational  time.  We  first  need  to  plot  the  mode  chart  along  with  the 
normalized  waveguide  frequency.  This  will  indicate  the  modes  present  in  the  structure. 
Several  iterations  were  done  to  obtain  a  viable  structure.  The  results  are  shown  in  Figure 
6.4.  This  corresponds  to  a  waveguide  with  the  following  parameters: 


B-V  Curve 


Figure  6.4  Mode  chart  for  the  single  mode  rectangular  waveguide  using  the 
Marcatili’s  model,  waveguide  width  2.0  pm 
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b  =  2  um  J=0.8  Lim  n\=  2.362  «2  =  «3  =  2.34 

,  (  6-6  ) 
«4  =  «5  =  1 

where  Figure  6.3  was  used.  Since,  in  Marcatili’s  method,  the  waveguides  are  assumed  far 
from  cutoff,  the  waveguide  structure  was  chosen  to  between  the  two  lower  modes,  Vwg.  A 
waveguide  with  a  width  of  1.5  pm,  which  is  closer  to  the  cut  off,  was  also  modeled  using 
the  Marcatili’s  method.  The  results  for  the  smaller  width  must  be  used  carefully,  since 
Marcatili’s  method  is  not  reliable  close  to  cut  off.  The  results  of  this  simulation  are 
shown  in  Figure  6.5.  For  the  height  of  the  cladding  region  (  n2  and  ),  it  was  desired  that 
the  evanescent  wave  would  be  sufficiently  attenuated  to  reduce  the  losses  due  to  surface 
roughness  at  the  clad-air  and  the  clad-sapphire  interfaces.  The  evanescent  waves 
attenuates  as: 

(6-7) 


B-V  Curve 


Figure  6.5  Mode  chart  for  the  single  mode  rectangular  waveguide  using  the 
Marcatili’s  model,  waveguide  width  1.5  pm 
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starting  at  the  core-clad  interface  and  into  regions  2  and  3.  In  the  equation  above,  x  would 
be  the  distance  from  the  interface  and  72,3  are  given  by  Equations  (  2-39  )  and  ( 2-40  ).  An 
additional  constrain  is  the  total  height  of  the  epitaxial  layer  that  can  be  grown  on  the 
sapphire  with  the  in-house  MBE  machine,  which  is  3  pm.  A  1.2  pm  layer  was  chosen  for 
the  cladding  in  contact  with  the  sapphire  and  1  pm  for  the  cladding  in  contact  with  air. 

The  cladding  in  contact  with  sapphire  was  chosen  to  be  larger  because  the  dislocation 
density  gets  reduced  as  you  move  away  from  the  sapphire/AlGaN  interface.  At  these 
thicknesses  the  attenuation  of  the  electric  field  is  (  using  Equations  (  2-39  ),  (  2-40  ),  and 
(  6-7  )  )  0.041  for  1  pm  and  0.022  for  1.2  pm.  From  this,  the  design  of  the  waveguide 
structure  is  given  in  Figure  6.6.  The  minimum  coupling  length  needed  for  complete 
power  transfer  for  the  2.0  pm  width  is,  using  Figure  6.3  and  Equations  (  2-41  ),  (  4-12  ) 
and  (  4-13  ),  ~  483  pm,  and  for  the  1.5  pm  width  the  coupling  length  is  -209  pm. 
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Figure  6.6  Side  view  of  the  single  mode  waveguide  structure  for  AlGaN,  x  is  chosen  to 
be  1 .5  and  2.0  pm 
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After  using  Marcatili’s  model  the  design  for  a  directional  coupler  with  2.0  pm 
waveguides  was  further  refined.  BeamProp  from  Rsoft  Inc  was  used  as  the  beam 
propagation  method  program.  Instead  of  creating  a  program  from  scratch  the  commercial 
program  was  used.  For  the  directional  coupler  design,  a  straight  input  was  used  to 
simplify  the  design  and  reduce  the  computational  time,  see  Figure  6.7.  This  enables  one 
to  concentrate  on  the  interaction  region.  In  Figure  6.7  one  sees  a  cross-sectional  top  view 
of  the  intensity  profile  of  the  optical  signal  in  the  directional  coupler.  The  cross-section 
was  taken  at  half  the  height  of  the  core  region.  In  it  one  can  see  the  light  coupling  back 
and  forth  between  the  guides.  Figure  6.8  is  the  power  in  the  core  for  each  waveguide  path 


Transverse  Direction  (pm) 


2.0 


0.0 


Figure  6.7  Intensity  profile  of  a  straight  directional  coupler,  Waveguide  width  2.0  pm  Based 
on  Brunner  et  al..1] 
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Power  transfer  between  the  waveguides 


Figure  6.8  Mode  power  for  each  waveguide,  waveguide  width  2.0pm.  Based  on  Brunner 
et  al..11 

as  a  function  of  the  propagation  direction.  This  figure  is  a  much  better  illustration  of  how 
the  optical  signal  couples  between  the  waveguides.  Figure  6.8  tends  to  suggest  the 
presence  of  an  additional  mode  in  the  structure.  Figure  6.9’s  cross  sectional  front  view  of 
the  directional  coupler  shows  a  weak  third  mode,  according  to  the  beam  propagation 
method,  between  the  guides.  This  third  mode  is  not  present  in  the  Marcatili’s  method 
analysis,  as  illustrated  in  Figure  6.10.  Furthermore,  Figure  6.8,  gives  the  minimum  length 
for  complete  power  transfer  between  the  guides  as  ~89  pm  which  is  much  smaller  than 
what  Marcatili’s  theory  predicts,  ~  483  pm. 
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Figure  6.9  Cross  sectional  view  of  the  directional  coupler  in  the  direction  of  propagation, 
according  to  the  beam  propagation  method  by  BeamProp.  Waveguide  width  is  2.0pm. 
Based  on  Brunner  et  al..71 
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B-V  Curve 


Figure  6.10  Mode  chart  for  the  interaction  region  of  the  directional  coupler,  using 
the  Marcatili’s  model,  corresponding  to  the  2  pm.  Total  width  is  4  pm.  Based  on 
Brunner  et  al.  71 
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Power  transfer  between  the  waveguides 


Figure  6.11  Mode  power  for  each  waveguide,  Waveguide  width  1.5pm.  Based  on 
Brunner  et  al .  ,71 

The  directional  coupler  with  the  1.5  pm  branches  was  simulated  with  the  beam 
propagation  method,  the  results  are  shown  in  Figure  6.1 1.  With  this  width  the  third  mode 


B-V  Curve 


Figure  6.12  Mode  chart  for  the  interaction  region  of  the  directional  coupler,  using  the 
Marcatili’s  model,  corresponding  to  the  1.5  pm.  Total  width  is  3  pm.  Based  on  Brunner 
et  al..11 
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Sapphire  orientation 
AIGaN/GaN  orientation 
-large  flat 


Small  Flat 


Figure  6.13  Crystal  orientation  between  sapphire  and  AlGaN 
appears  to  disappear.  There  are  still  discrepancies  between  the  simulation  methods.  In  the 
Marcatili’s  method  there  is  no  third  mode  and,  according  to  Figure  6.12,  the  second  mode 
is  close  to  cut  off.  Additionally,  it  gives  a  coupling  length  for  complete  power  transfer  of 
~209  pm  while  the  beam  propagation  method  gives  us  ~  53  pm.  From  Figure  6.1 1  we 
can  see  that  the  third  mode  was  not  entirely  gone  in  the  first  1000  pms,  it  dissipated 
slowly. 

First  Mask  Set  Processing,  and  Results 

With  these  analyses  two  initial  mask  templates  were  designed.  In  it  there  were 
waveguides  of  different  widths  and  lengths  along  with  different  configurations  of 
directional  couplers,  see  APPENDIX  A.  These  mask  templates  were  done  to  produce 
solely  passive  devices,  initiate  the  development  of  device  processing  in  AlGaN,  and 
determine  which  structures  are  suitable  for  a  directional  coupler  design.  The  difference 


75 


between  the  two  mask  templates  was  that  on  one  of  them  the  structures  were  defined  by 
the  metal  in  the  template  (light  field  mask)  and  the  other  by  the  absence  of  it  (  dark  field 
mask  ).  The  AlGaN  material  was  grown  in-house  on  an  MBE  machine  using  sapphire 
substrates.  One  disadvantage  of  using  the  sapphire  substrates  is  that  the  crystal  orientation 
of  sapphire  is  30°  offset  to  that  of  the  AlGaN  material  system,  see  Figure  6.13.  This 
makes  it  very  difficult,  if  not  impossible,  to  cleave  the  sample.  A  good  method  to  cut  the 
sample  is  to  saw  it  all  the  way. 

Being  a  material  that  withstands  high  temperature  and  caustic  environments  makes  it 
difficult  to  etch.  Up  to  date  the  most  successful  methods  to  etch  AlGaN  are  dry  etching, 
although  some  progress  has  been  done  with  wet  etching  techniques.  ’’’’’’  For  this 
work,  the  load  locked  Plasma-Therm  SFR  770  ICP  reactor  was  used  to  etch  AlGaN. 
Initially,  there  was  not  a  suitable  mask  or  recipe  to  etch  a  waveguide  that  was  deep 
enough  or  steep  enough.  For  this  particular  application  the  waveguides  had  to  be  ~  2  pm 
high  with  steep  sidewalls,  see  Figure  6.6.  Many  combinations  of  masks  and  recipes  were 
investigated  until  a  suitable  one  was  found,  see  APPENDIX  B.l.  A  suitable  masking 
material  turns  out  to  be  Ni.  This  was  the  only  material  able  to  withstand  the  etching  of  the 
AlGaN  materials  and  have  slow  degradation.  The  other  materials  had  to  be  at  least  half 
the  thickness  of  the  desired  etching  height.  Standard  processes  were  used  as  the  initial 
starting  points  to  transfer  the  pattern  into  the  AlGaN  material.  The  lift  off  process  was 
used  to  create  a  metal  pattern  while  the  etch  masking  process  was  use  for  a  photoresist 

'l  r 

pattern  with  the  dark  and  light  field  masks  respectively.  Figure  6.14  shows  the 
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Light  Field  Mask  Template  Processing 


Dark  Field  Mask  Template  Processing 
(Lift-off  Processing) 


Evaporate  metal 
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finished  for  the  dark 
field  processing. 
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PMGI-SF  11  / 

AlGaN 

PMGI-SF  11 

AlGaN 

Figure  6. 14  Process  flow  for  the  transfer  of  the  Mask  template  pattern  into  the  sample 


schematic  flow  of  the  process  followed  for  both  the  dark  field  and  the  light  field 
templates.  APPENDIX  C.l  shows  the  details  of  the  dark  field  process.  The  process  in 
APPENDIX  C.l  adds  an  optional  SiN  layer  to  detach  Ni  at  the  end  of  the  process.  In  both 
processes,  a  dual  layer  of  positive  photoresist  is  used  (1813  and  SF-1 1). 

Although  sawing  the  sample  is  effective  in  dividing  it  into  smaller  pieces  it  doesn’t 
leave  the  end  points  of  the  waveguides  in  a  good  enough  condition  to  launch  light  into  or 
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Figure  6.15  End  polishing  of  a  waveguide  using  the  FIB 


out  of  them.  To  polish  the  ends  of  the  waveguides  a  FEI  Dual  Beam  620  FIB/SEM 
workstation  was  used.  This  type  of  polishing  has  been  shown  to  improve  the  facets  of 
GaN  laser  diodes.81  Figure  6.15  shows  the  results  of  polishing  with  the  Focused  Ion  Beam 
(FIB).  As  it  can  be  seen,  the  polishing  of  the  end  of  the  waveguide  is  much  better  than  the 
etching  from  the  ICP. 

Another  problem  that  showed  up  was  that  once  the  light  was  launched  into  the 
waveguide  the  output  was  obscured  by  the  background  stray  light  emanating  from  the 
laser.  In  order  to  solve  this,  a  -45°  angle  was  cut  into  the  waveguide  with  the  FIB,  this 


(a)  (b) 


Figure  6.16  Etched  mirror  (a)  close  up  of  the  mirror,  (b)  mirror  with  respect  to  spec 
present  near  the  waveguide 
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Figure  6.17  Schematic  of  testing  setup 

made  the  light  shoot  straight  up  thus  avoiding  any  stray  light  from  the  laser  (  see  Figure 
6.16  ).  The  mirror  shown  on  Figure  6.16  was  etched  at  ~  1mm  from  the  input  end  of  the 
waveguide. 

The  schematic  of  the  experimental  setup  used  to  launch  light  into  the  waveguides  and 
examine  them  from  the  top  is  shown  on  Figure  6.17.  In  the  experimental  setup  a  5  mW, 
632.8  nm,  laser  light  is  expanded  and  then  collimated,  not  shown.  The  laser  light  is  then 
passed  through  a  polarizer  in  which  only  light  polarized  normal  to  the  surface  of  the 
sample  is  selected,  see  Equation  (  6-5  ).  The  light  is  then  focused  into  the  waveguide. 
Alignment  is  made  through  a  5  degrees  of  freedom  stepper  motor  system.  The  stepper 
motor  system  can  be  controlled  by  a  computer  or  manually.  Light  is  then  collected  from 
the  top  of  the  sample  with  a  microscope  system  and  imaged  into  a  camera.  The  camera 
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(a)  (b) 


Figure  6.18  Light  guidance.  Both  images  are  of  the  same  region  (a)  output  with  lights  on 
and  (b)  output  with  lights  off 

was  a  COHU  6310  series  monochrome  camera  capable  of  640  x  480  resolution  and  the 
output  is  fed  into  an  image  capture  card. 

Figure  6.18  shows  the  output  from  the  waveguide  shown  in  Figure  6.16.  This  particular 
waveguide  was  ~6  pm  wide,  and  thus  multimode  on  both  the  simulation  models,  but  it 
illustrated  the  success  of  the  wave  guiding  capability  of  the  structure  and  the  success  of 
the  ~45°  mirror.  Light  escaping  from  the  top  of  the  waveguide  (not  the  mirror  output) 
cannot  be  observed,  although  a  very  strong  signal  is  emerging  from  the  waveguide.  Figure 


Figure  6.19  Light  coupling  into  the  spec  shown  in  Figure  6.16  b) 
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6.19  shows  the  spec  next  to  the  waveguide  lighted  up,  see  Figure  6.16  b).  It  is  suspected 
that  most  of  the  loss  mechanism  in  the  waveguide  is  through  the  rough  sidewalls  where 

the  core  is  exposed  to  air,  see  Figure  6.15  and  Figure  6.16.  It  is  known  that  interface 

22 

deformations  will  result  in  losses. 

In  order  to  see  if  the  light  was  escaping  the  waveguides  from  the  sides  rather  than  from 
the  top,  a  modification  was  made  to  the  setup  to  see  a  greater  portion  of  the  sides.  The 
sides  were  viewed  at  an  angle  of  ~  11°  from  normal  and  the  amount  of  light  observed  was 


c)  d) 


Figure  6.20  Light  being  launch  in  a  directional  coupler  and  viewed  a),  c)  at  11°  from 
normal,  b),  d)  viewed  normal  to  the  sample,  c)  and  d)  were  taken  1 .2  mm  down  from  a) 
andb) 
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a)  b) 

Figure  6.21  Output  of  the  directional  coupler,  a)  electron  micrograph  b)  optical  output 
greater  than  looking  at  normal  incidence.  The  device  being  viewed  from  the  sides  and  top 
was  a  passive  zero-gap  directional  coupler,  see  Figure  6.20,  with  the  interaction  and 
branches  regions  having  a  width  of  ~  5  pm  and  ~  3  pm  respectively.  The  amount  of  light 
observed  was  still  dim  and  was  difficult  to  analyze.  At  the  output  both  branches  had  an 
optical  signal  present  at  a  distance  of  ~  2  mm  from  the  front  end,  Figure  6.21.  The  signal 
had  equal  intensity.  Due  to  the  widths  of  the  waveguides,  they  were  multimode. 

A  problem  that  occurs  is  that  the  widths  of  the  devices  using  the  lift-off  process,  see 
Figure  6.14,  tends  to  be  bigger  than  the  structures  in  the  mask  template.  In  order  to 
compensate  for  this,  a  structure  with  5%  A1  concentration  was  grown,  see  Figure  6.6. 
New  devices  were  grown  with  this  structure.  Figure  6.22  shows  the  result  of  the  device 
constructed,  a  bright  signal  light  is  shown  coming  from  the  end  of  the  waveguide.  The 
problem  is  that,  according  to  calculations  (  see  Figure  6.23  )  the  waveguide  should  not 


have  been  guiding  at  all. 
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Figure  6.22  1.4  pm  waveguides  using  AlGaN  at  5%  A1  concentration  for  cladding  layer 
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Figure  6.23  Simulation  of  a  1.4  pm  waveguide  by  the  beam  propagation  method. 
Based  on  Brunner  et  alJx 
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Figure  6.24  Indices  of  refraction  based  on  Bergman  et  al.  at  room  temperature  and  632.8 
nm,  reference  69 

In  trying  to  obtain  better  agreement  with  the  results,  another  index  of  refraction  model 
for  AlGaN  was  used.69  This  model  was  obtained  by  M.  J.  Bergmann  et  al.,  see  Chapter  5. 
From  equations  (  5-2  )  -  (  5-4  )  and  Table  2  we  obtain  the  refraction  index  graph 
illustrated  on  Figure  6.24.  For  the  Al  concentration  of  interest  we  get,  at  room 
temperature: 


n0  ne 
2.268  2.304 
2.236  2.270 
2.254  2.289 


Core,  GaN 

Cladding,  Alo.1Gao.9N 
Cladding,  Al0.05Ga0.95N 
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The  new  refraction  indices  increase  the  index  difference  between  the  core  and  the 
cladding.  This  index  difference  is  given  by: 

^  core  ^cladding 

-  ( 6-8 ) 

ncox  e 

This  equation  gives  us  an  index  difference  of  0.0042  for  the  Brunner  et  a/.parameters71 
and  0.0065  for  the  Bergmann  et  al.69  parameters.  The  increase  index  difference  means 
that  the  modes  will  be  more  strongly  guiding.  These  were  enough  to  have  the  beam 
propagation  simulations  of  the  1.4  pm  waveguide  show  light  confinement  and  guidance, 
as  Figure  6.25  shows.  The  Marcatili  simulations  show  the  waveguide  in  cut  off,  non¬ 
guiding  (  see  Figure  6.26  ). 


2.0 


0.0 


X  (/am) 


Pathway, 

Monitor: 

1 ,  WG  power 


Figure  6.25  New  simulations  of  the  waveguide  that  shows  strong  guidance.  The 
parameters  for  this  one  is  as  follows:  width  1.4pm,  core  0.8pm,  Al  concentration  5%. 
Based  on  Bergmann  et  al.69 
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B-V  Curve 


Figure  6.26  Marcatili’s  simulation  using  the  new  refraction  indices  and  the  1.4  gm 
waveguide,  A1  concentration  5%.  Based  on  Bergmann  et  al..69 


B-V  Curve 


Figure  6.27  Marcatili’s  simulation  using  the  new  refraction  indices  and  the  2.0  gm 
waveguide,  Al  concentration  5%.  Based  on  Bergmann  et  al..69 
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More  precise  measurements  of  the  waveguide’s  dimension  were  made  using  an  Atomic 
Force  Microscope  (  AFM  ).  The  measured  width  of  the  waveguide  was  ~  2pm.  Also  the 
A1  concentration  is  not  exact,  according  to  the  MBE  grower  it  might  be  off  by  1-  2%. 
With  these  more  precise  measurements  more  simulations  were  made.  The  simulations 
showed  wave  guiding.  Figure  6.27  shows  one  possible  scenario  for  wave  guiding  using 
the  Marcatili  method.  All  of  this  shows  that  the  index  of  refraction  formulas  established 
by  M.  J.  Bergmann  et  al.  are  more  suited  for  the  design  of  the  waveguides,  at  least  for  the 
current  in-house  growth  conditions. 

Further  simulations  were  done  that  concentrated  on  the  coupling  region  Figure  6.28 
shows  the  Marcatili  simulation  for  the  coupling  region  of  a  directional  coupler  whose 
waveguides  are  2.0  pm  each,  4.0  pm  total  width  in  the  coupling  region.  According  to  the 


B-V  Curve 


Figure  6.28  Mode  chart  for  the  interaction  region  of  the  directional  coupler,  using  the 
Marcatili’s  model,  corresponding  to  the  2  pm.  Total  width  is  4  pm.  Based  on  Bergmann 
et  al.  ,69 
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simulation  there  are  two  modes  present  with  one  close  to  cut  off.  According  to  Marcatili, 
the  coupling  region  in  the  1.4  pm  directional  coupler  has  only  one  mode  present.  Also, 
the  Marcatili  simulations  gives  us  the  minimum  coupling  length  for  complete  power 
transfer  for  a  1.4  pm  waveguide  as  ~1  63  pm  and  for  the  2.0  pm  waveguide  the  coupling 
length  is  -458.4  pm.  With  the  new  index  of  refraction  the  attenuation  given  by  Marcatili 
is  ~  0.074  at  a  1pm  distance. 

For  the  beam  propagation  method  an  entire  directional  coupler  was  simulated.  Its 
results  are  shown  in  Figure  6.29.  An  advantage  that  the  beam  propagation  method  has 
over  the  Marcatili  method  is  that  of  being  able  to  simulate  the  entire  directional  coupler 
structure.  The  beam  propagation  method  shows  a  coupling  length  for  complete  power 


Pathway, 
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Figure  6.29  An  entire  directional  coupler  simulated  by  the  beam  propagation  method.  The 
waveguide  branches  are  1.4  pm  wide.  The  coupling  region  is  3750  pm  long. 
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transfer  of  only  63  pm  for  the  1.4  waveguide  width.  This  is  in  contrast  with  the  ~  163  pm 
with  the  Marcatili  simulation.  From  Figure  6.29,  note  that  once  in  the  coupling  region, 
there  is  a  transition  period  before  the  modes  settle  into  a  nice  two-mode  transition. 

Through  out  the  course  of  this  research,  small  variations  were  made  to  the  process  in  an 
attempt  to  improve  the  sidewall  roughness.  A  variation  to  the  ICP  etching  recipe  was  also 
attempted  with  better  results.  APPENDIX  B.2  shows  the  experimental  attempts  to 
improve  the  sidewalls.  The  logic  of  the  attempt  is  as  follows:  In  the  ICP  etching  there  is  a 
chemical  component  and  sputtering  component.  The  physical  sputtering  is  dominated 
by  the  accelerations  of  energetic  ions  in  the  plasma  towards  the  sample  material.  In 
reducing  the  rf  power,  the  contribution  of  the  chemical  etching  is  increased  while  that  of 
sputtering  is  reduced.  With  the  experimental  procedure  it  was  found  that  a  reduction  of 
the  rf  power  to  75  W  improved  the  sidewall  roughness,  with  little  damage  to  the  sample 
and  little  effect  to  the  sidewall  slope. 

During  the  experiments,  it  was  found  that  good  waveguide  widths  on  the  mask  template 
were  1.4  pm  and  2.0  pm.  This  is  due  to  processing  which  will  transfer  to  a  pattern  with 
larger  line  sizes.  Additionally,  the  separation  chosen  at  the  input  and  output  of  the 
waveguide  branches  was  too  small.  Figure  6.20  a)  and  b)  shows  that  when  light  was 
attempted  to  be  launched  into  one  of  the  input  branches  of  the  directional  coupler,  the 
spot  size  was  large  enough  to  couple  light  into  the  other  input  branch.  This  would  make  it 
difficult  to  evaluate  the  directional  coupler,  since  there  could  be  two  simultaneous  inputs 
to  the  coupling  region,  see  Figure  4.6. 
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Second  Mask  Set  Design,  Processing,  and  Results 
Based  on  this  information,  a  new  directional  coupler  was  designed,  see  Figure  6.30. 
This  new  design  for  the  directional  coupler  will  be  for  widths  of  1.4  and  2.0  pm.  In  this 
new  design  the  waveguides  in  the  input  sections  are  separated  by  100  pm  while  in  the  exit 
section  they  are  separated  by  40  pm.  This  is  done  so  that  the  output  sections  are  not  in  the 
direct  path  of  the  input  optical  beam.  The  bending  sections  were  designed  to  have  an 
angle  of  0  ~  6°  (see  Figure  4.6)  although  Equation  (  4-32  )  will  give  us  5.6°  and  3.9°  for 
the  1.4  and  2.0  pm  widths  respectively.  The  difference  in  the  entrance  angle  stems  from 
the  fact  that  it  is  the  instantaneous  angle  that  needs  to  meet  Equation  (  4-32  )  criteria  and 
not  the  total  angle,  see  Chapter  IV.  For  the  second  mask  set,  the  coupling  region  was 


500  pm 


Figure  6.30  New  directional  coupler  design. 
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Ohmic  Contacts 


AlGaN 
GaN 

AlGaN 

Figure  6.31  Metal  contacts  on  the  directional  coupler, 
varied  from  3000  to  3200  pm  in  increments  of  10  pm,  while  maintaining  the  same  total 
length  of  the  device.  This  would  enable  us  to  determine  the  experimental  coupling  length 
and  compare  it  with  the  models. 

For  the  second  mask  template,  that  encompassed  the  structure  on  Figure  6.30, 
additional  attempts  were  made  to  create  an  active  device.  As  shown  in  Figure  6.1,  the 
structure  to  be  envisioned  would  be  an  n-(semi-i)-n  and  if  possible  an  p-(semi-i)-n. 
Contacts  were  placed  as  shown  in  Figure  6.31.  The  bottom  layer  served  as  a  common 
ground  for  the  devices.  This  will  enable  the  electric  field  to  flow  from  one  contact  to  the 
other.  By  having  the  AlGaN  conductive  and  the  GaN  as  insulating  as  possible  it  was 
hoped  that  the  electric  field  would  be  concentrated  in  the  GaN  core  region.  In  creating  a 
process  to  accomplish  this,  several  factors  have  to  be  taken  into  account.  One  of  them  is 
the  placement  of  the  metal  on  top  of  the  coupling  region.  The  scope  of  this  problem 
entails  the  precise  placement  of  a  metal  line  on  top  of  a  mesa  -3000  pm  long  and  -2.8 
pm  wide  with  little  room  to  spare.  This  needed  to  be  done  consistently  through  out  the 
wafer.  Prior  to  the  complete  mask  development,  tests  were  made  to  rule  out  alternatives. 
One  attempt  was  placing  the  top  metal  contact  before  the  AlGaN  etching.  This  was  done 
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Figure  6.32  Test  for  placing  a  metal  contact  on  top  of  the  coupling  region  in  a  directional 
coupler. 

with  the  hope  that,  during  the  etching,  any  metal  contact  not  aligned  with  the  coupling 
region  would  be  etched  out.  The  problem  was  that  Ni  mask  and  a  possible  SiN  under 
layer  were  removed  with  chemicals  that  would  attack  the  ohmic  contacts.  The  ohmic 
contacts  used  for  the  n-type  region  were  Ti/Al/Ti/Au,  with  the  Ti  in  contact  with  the 
surface.  After  depositing  the  metals,  the  sample  goes  through  a  rapid  thermal  annealing 
(  RTA  )  cycle  to  create  an  alloy.  The  Ni  mask  could  not  be  left  on  the  structure  because  it 
would  apply  an  electric  field  through  out  the  structure  and  it  could  be  an  avenue  for  other 
problems  to  emerge.  The  one  that  was  successful  was  the  placement  of  the  metal  contact 
after  the  directional  coupler  mesa  was  etched,  as  shown  in  Figure  6.32. 

Another  problem  was  the  placement  and  connection  of  a  pad  that  would  make  contact 
to  the  metal  on  top  of  the  mesa.  Due  to  the  dimensions  of  the  coupling  region  it  was 
impossible  to  place  an  electrical  probe  on  it.  A  separate  metal  pad  had  to  be  created  that 
would  make  contact  with  it.  Implantation  of  an  insulating  region,  or  any  type  of 
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Figure  6.33  Ohmic  contact  configuration  with  metal  pad  and  an  air  bridge 


implantation  in  AlGaN  for  that  matter,  was  not  an  option.  The  approach  taken  was  to  etch 
a  trench  all  the  way  to  the  sapphire,  which  is  insulating  in  nature.  An  air  bridge  was  used 
to  make  contact  between  the  metal  pad  and  the  mesa  contact.  Figure  6.33  illustrates  this 
approach.  The  metal  used  for  the  bridge  and  the  pad  at  the  bottom  of  the  trench  didn’t 
have  a  special  contact  or  electrical  requirement  to  be  met. 

With  this  information  a  second  mask  template  set  was  designed.  This  one  was  targeted 
to  further  investigate  the  passive  directional  coupler  as  well  as  to  develop  the  processing 
steps  necessary  to  create  an  active  directional  coupler.  The  steps  that  were  to  be 
implemented  are  shown  in  Figure  6.34.  Figure  6.34  shows  the  crossectional  view  as  seen 
from  the  coupling  region.  The  ohmic  contacts  to  be  deposited  on  top  of  the  directional 
coupler  and  on  the  common  ground  plane  are  identical  for  the  n-i-n  structure.  But,  on  a  p- 
i-n  structure  they  would  be  different,  therefore  different  mask  templates  had  to  be  used 
for  each  one.  For  the  p-type  contact,  the  most  often  used  alloy  is  Ni/Au,  which  was  run 
through  an  RTA  cycle  to  obtain  the  desired  ohmic  contact. 
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Figure  6.34  Overview  of  the  processing  steps  that  were  implemented  on  the  second 
mask  set. 
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Figure  6.35  Top  view  diagram  of  the  active  directional  coupler.  Not  to  scale 
Figure  6.35  shows  the  top  view  of  the  active  directional  coupler  building  cell.  The 
interleaving  fingers  are  the  bridges  shown  in  Figure  6.33.  Instead  of  using  a  single  bridge 
multiple  smaller  bridges  were  used.  This  made  it  easier  to  take  out  the  photoresist  from 
underneath  the  bridges.  Following  the  processing  steps  of  Figure  6.34,  each  of  the 
different  shades  indicates  a  separate  feature  that  was  transferred  from  a  mask  template. 
There  were  a  total  of  six  mask  templates  for  the  second  set.  Of  the  six  five  are  shown  in 
Figure  6.35.  Additionally,  redundancy  was  built  in  the  mask  set  by  having  sets  of  four 
directional  couplers  with  equal  coupling  length,  before  increasing  to  the  next  length. 

With  the  new  mask  set  some  problems  were  more  prevalent  in  constructing  the 
directional  coupler.  One  of  the  problems  was  at  the  waveguide  junction  in  the  directional 
coupler,  Figure  6.36  a).  As  the  waveguides  come  closer  together  the  width  of  the  material 
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a) 


b) 

Figure  6.36  Waveguide  Junction  in  a  directional  coupler  a)  sketch  of  the  template  mask 
for  the  1.4pm,  b)  failed  transfer  of  the  pattern  to  the  photoresist. 

separating  them  becomes  smaller  and  smaller.  In  theory  it  becomes  very  small  until  it 
vanishes.  In  actual  manufacturing  the  minimum  separation  is  dependent  on  the  mask 
template  writing  tools  and  lithography  resolution.  For  the  second  mask  set,  the  resolution 
was  -0.25  pm.  For  the  transferred  pattern  the  process  itself  will  also  weigh  in  on  the  final 
feature  size. 

For  the  pattern  transfer  using  the  lift-off  processing,  see  Figure  6.14,  dual  layer  photo 
resist  was  used.  When  attempting  to  establish  the  pattern  into  the  photo  resist  the  material 
between  the  waveguides  peeled  off  as  shown  in  Figure  6.36  b).  This  presents  a  problem 
when  constructing  the  waveguide,  since  it  will  create  an  opening  in  the  path  of  the  optical 
signal  and  disrupt  it.  The  source  of  this  problem  can  be  seen  in  Figure  6.37.  As  seen  in 
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Figure  6.37  Crossectional  view  of  the  photoresist  pattern  transfer  of  the  junction  region 
of  the  directional  coupler. 

Figure  6.14  and  Figure  6.37  for  the  metal  lift-off,  the  bottom  photoresist  layer  etches 
slightly  more  than  the  top  layer.  This  is  convenient  since  otherwise  the  metal  would  have 
been  deposited  along  the  sidewall,  in  addition  to  the  bottom  of  the  pattern.  This  would 
create  winglets  in  the  metal  pattern,  which  would  interfere  with  the  pattern  transfer  to  the 
AlGaN/GaN  material.  This  works  against  the  shrinking  separation  at  the  junction. 
Looking  at  Figure  6.37  we  can  see  that  at  a  certain  point  the  bottom  photoresist  layer 
would  detach  itself  from  the  AlGaN  material  while  leaving  the  top  photoresist  layer,  thus 
the  peeling  as  seen  in  Figure  6.36  would  occur. 

When  a  mask  template  is  produced  the  curves  are  not  created  in  a  continuous  fashion. 
Rather,  the  curves  are  created  in  a  stepwise  fashion,  as  is  seen  in  Figure  6.36.  This  can  be 
used  to  our  advantage  in  creating  a  process.  In  order  to  minimize  the  peeling,  the  process 
has  to  be  tweaked  so  that  the  photoresist  between  the  waveguides  is  not  etched  enough  to 
peel  or  etch  away  the  entire  finer  step  and  leave  the  next  step  increase  intact.  To  do  this 
the  exposure  times,  the  developing  times,  and  the  O2  plasma  asher  times  had  to  be 
tweaked  to  obtain  the  desired  characteristics.  The  process  follower  in  APPENDIX  C.2 
for  the  “WG  definition”  process  shows  the  result  of  these  iterations. 
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Figure  6.38  SEM  of  etched  samples  with  uneven  etching 


With  the  new  template  mask  set,  a  few  samples  were  etched.  When  the  samples  were 
etched  they  had  a  dull  black  appearance  to  the  naked  eye.  When  a  Scanning  Electron 
Microscope  (SEM)  was  used  it  was  evident  that  the  samples  were  etched  unevenly, 
unetched  pillars  were  permeating  the  entire  exposed  region,  see  Figure  6.38.  Some  of  the 
previous  samples  turned  out  this  bad  and  some  to  a  lesser  degree.  This  behavior  was 
investigated  to  determine  its  causes.  A  small  design  of  experiments  (  DoE  )  was 
conducted  towards  this  end,  see  Table  4.  These  studies  showed  that  there  is  an  interaction 
between  the  freon  etch  (F-23  or  FI 4)  that  seeds  the  formation  of  the  column  structures  in 
the  ICP  etch.  If  the  AlGaN/GaN  sample  with  no  S13N4  was  put  through  the  freon  process, 
the  ICP  etch  comes  out  clean  (no  columns).  This  is  also  true  if  a  sample  with  S13N4  is 


Table  4  Small  DOE 
SiN 

thickness 

in 


Trials 

Sample 

Freon  (RIE) 

SiN 

angstroms 

Asher 

Pattern 

mask 

Output 

1 

252-2-1 

Over  etched 

no 

no 

no 

none 

Good 

2 

252-2-2 

Over  etched 

yes 

5751 

no 

no 

none 

Bad 

3 

252-2-3 

no 

no 

yes 

no 

none 

Good 

4 

252-2-4 

Over  etched 

yes 

5751 

yes 

no 

none 

Bad 

5 

252-2-5 

Over  etched 

yes 

5751 

yes 

yes 

SiN 

Bad 
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etched  without  the  freon  step.  The  freon  step  is  used  to  pattern  the  ShN4  for  later  etching 
the  AlGaN/GaN  sample.  Over  etching  of  the  ShN4  was  done  with  little  or  no  success,  the 
freon  etch  doesn’t  etch  the  GaN/AlGaN  and  the  surface  is  still  smooth  at  the  end.  When 
the  sample  was  over  etched  an  ED  AX  (Energy  Depressive  X-ray  Spectrum)  analysis  was 
done  (see  Figure  6.39).  This  revealed  that  there  was  fluorine  residue  in  the  sample.  The 
analysis  only  shows  elements  presents  and  not  chemical  formulas. 

It  is  seems  that  the  fluorine  residue  by  itself  or  in  a  chemical  composition  with  another 
element  present  prevents  the  etching  of  the  AlGaN/GaN  material.  Several  attempts  in 
cleaning  the  sample  after  a  freon  etch,  but  before  the  ICP  etch,  were  done.  This  showed 
that  by  over  etching  the  ShN4  and  then  dipping  the  sample  for  5  min  in  phosphoric  acid  a 
clean  ICP  etch  was  obtained.  This  step  does  not  work  if  there  is  ShN4  remaining,  before 
the  rinse,  on  the  area  intended  for  etching.  BOE  could  also  be  used,  but  since  it  dissolves 
ShN4  it  needs  to  be  diluted  and  have  a  controlled  timing  to  etch  the  affected  area. 
Additionally,  it  etches  ShN4  isotropically  and  it  would  under  cut  the  patterns.  This 
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Figure  6.39  ED  AX  of  sample  after  freon  over  etch 
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Figure  6.40  Directional  coupler  output  region 
prompted  a  slight  change  in  the  process.  RIE  cleaning  of  SfdSU  was  eliminated.  Now 
either  no  Si3N4  was  used  or  it  was  left  to  be  patterned  with  the  AlGaN  in  the  ICP,  the 
optional  parts  in  the  “WG  definition”  process  follower  in  APPENDIX  C.2. 

Once  the  process  was  de-bugged,  new  devices  were  constructed  using  the  5%  A1 
concentration.  Figure  6.40  shows  the  output  region  of  the  directional  coupler.  This  figure 
clearly  shows  the  angle  cut  on  the  waveguides  plus  an  angle  cut  between  the  waveguides. 
This  angle  cut  between  the  waveguide  was  done  in  order  to  test  for  any  wave  guiding  that 
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Figure  6.41  Launching  convention  for  the  experiments 


a)  Left  side  launch  b)  Right  side  launch 

Figure  6.42  Light  exiting  the  directional  coupler  with  a  theoretical  interaction  length  of 
3190  pm.  Magnification  is  200x 

might  have  occurred  in  the  AlGaN  layer  directly  on  top  of  the  sapphire,  see  Figure  6.6. 
Additionally,  the  center  cut  helps  in  locating  the  output  section  using  the  testing  set  up  of 
Figure  6.17. 

Light  was  launched  into  the  left  side  and  then  into  the  right  side  of  the  device.  Figure 
6.41  shows  the  launching  conventions  for  the  experiments.  The  output  was  monitored 
using  the  camera  set  up  in  Figure  6.17.  Figure  6.42  shows  the  output  of  one  such  device 
in  which  light  was  launched  into  the  a)  left  side  b)  and  the  right  side  of  the  device.  In  both 
images  we  see  that  one  of  the  branches  is  brighter  than  the  other  one,  and  that  they  change 
when  changing  the  input  branch.  Theoretically  this  is  expected  due  to  the  symmetry  of  the 
device.  It  is  an  indication  that  the  device  might  be  working  properly,  at  least  in  the 
passive  mode.  This  simple  test  reduces  the  probability  that  at  least  one  of  the  branches 
might  be  obstructed  or  damaged.  This  might  present  a  problem  in  analyzing  the  structures 


101 


a)  b) 


Figure  6.43  Higher  magnification  (  800x  )of  the  directional  coupler  of  Figure  6.42  with  a) 
the  lights  off  and  b)  the  lights  on. 

because  it  will  skew  the  results.  Several  devices  of  different  lengths  were  tested  but  gave 
indications  of  partial  obstructions  in  the  waveguides. 

Additional  testing  of  the  waveguides  was  preformed.  Figure  6.43  shows  the  same 
directional  coupler  as  Figure  6.42  at  a  higher  magnification.  Figure  6.43  b)  shows  all 
three  angle-cuts  shown  in  Figure  6.40.  At  this  magnification  the  signals  were  so  strong 
that  they  appeared  to  be  undistinguishable  in  magnitude,  they  saturated  the  detector  (  see 
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Figure  6.44  Attenuation  of  the  laser  using  the  directional  coupler  of  Figure  6.43.  Light 
was  launched  on  the  left  side.  Light  launched  at  the  left  side. 
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Figure  6.45  Right  side  launching  in  same  structure  as  in  Figure  6.44  and  with  attenuation 

Figure  6.43  a) ).  Because  of  this,  an  attenuator  had  to  be  used  with  the  laser.  It  was  only 
then  that  the  intensity  difference  was  measurable,  see  Figure  6.44.  Figure  6.44  and  Figure 
6.45  are  for  left  and  right  launched  light  respectively.  The  measurements  were  obtained 
from  the  image  by  reading  the  intensity  of  the  pixel,  which  was  scaled  from  0-255,  for 
an  8  bit  gray  scale. 

In  Figure  6.44,  the  dual  peak  could  be  caused  by  an  imperfection  in  the  facet  milling, 
since  only  one  waveguide  showed  it,  and  the  branch  is  not  in  a  dual  mode  region.  For  the 
Marcatili  approximation,  the  results  are  skewed  because  we  have  a  tapered  entrance  and 
exit  from  the  device,  see  the  section  on  the  “bending  section”  in  Chapter  4.  Figure  6.44 
and  Figure  6.45  illustrates  the  problem  with  uneven  attenuation  in  the  output  branches. 
This  uneven  attenuation  could  have  been  caused  by  the  process,  resulting  in  small 
imperfections  in  the  waveguides.  This  would  render  many  of  the  waveguides  unusable  for 
a  direct  comparison.  In  order  to  circumvent  this  problem  some  mathematical 
manipulation  can  be  accomplished  in  order  to  extract  the  data  from  the  experiment. 
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If  we  assume  that  in  the  coupling  region  any  imperfection  or  attenuation  affects  both 
modes  equally,  or  an  insignificantly  small  uneven  amount,  then  the  output  branches 
would  be  the  source  of  uneven  attenuation.  Therefore,  if  the  output  ratio  of  the  intensity 
with  light  incident  on  one  of  the  branches,  say  the  right  one,  is  given  by: 


C.bP, 


aPx 


(6-9) 


where  a  and  b  are  the  attenuation  in  each  of  the  branches,  P\  is  the  ideal  output  power  in 
the  incident  branch  and  P2  is  the  ideal  output  power  in  the  non-incident  branch  (  without 
any  attenuation).  Then  output  power  ratio  with  the  incident  light  in  the  other  branch,  say 
the  left  one,  is  given  by: 

aP'-, 


D  = 


bp; 


(6-10) 


where  the  prime  indicates  light  incident  on  the  branch  opposite  to  Equation  (  6-9  ).  The 
attenuations  a  and  b  change  since  they  are  fixed  to  the  individual  branches.  By  the 
symmetry  of  the  device  we  have  Px  =  P{  and  P2  =  P' .  By  dividing  Equation  (  6-9  )  by 
Equation  (  6-10  )  we  obtain 


=  =  (6-11) 
D  a 2 

where  AR  would  be  the  attenuation  ratio.  With  this  value  we  can  now  obtain  the  power 
ratio  without  the  uneven  attenuation. 

A  full  set  of  passive  waveguides  were  grown  and  tested  with  the  experimental  set  up  of 
Figure  6.17.  In  testing  the  waveguides  there  was  some  vibration  present  in  the  set  up, 
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although  the  entire  set  up  was  mounted  in  vibration  isolation  tables.  For  testing  the 
waveguides  light  was  launched  in  either  the  left  side  or  the  right  side  of  the  input 
branches.  In  the  output,  several  images  were  taken.  Since  the  vibration  smudged  many  of 
the  images  only  those  that  were  not  smudge  were  selected.  Once  one  of  the  branches  was 
finished,  we  proceeded  to  the  next  one.  For  each  of  the  branches  several  images  were 
taken.  They  were  analyzed  in  the  same  fashion  as  in  Figure  6.44  and  Figure  6.45  and  the 
ratio  of  both  output  branches  was  obtained. 

All  the  data  concerning  a  branch  was  statistically  analyzed  by  averaging  and  obtaining  a 
standard  deviation.  Table  5  shows  a  sample  of  the  data  taken  for  a  single  directional 
coupler.  The  R  and  L  notation  signifies  the  output  power  from  the  right  or  left  branch 
respectively.  Using  Equation  (6-11  )  the  uneven  attenuation  was  then  canceled  out  of  the 
data.  Figure  6.46  shows  the  data  obtained  from  various  directional  couplers  of  the  same 
sample  and  the  corrected  values  of  the  data.  There  was  no  specific  pattern  in  the  uneven 
attenuation  implying  that  its  origin  was  not  in  the  material,  but  the  structure  itself.  In  two 
of  the  structures  only  a  single  branch  was  available  to  launch  the  light  into  and  therefore 
impossible  to  correct. 


Table  5  Sample  of  data  taken  from  a  single  directional  coupler.  The  ratio  is  set  as 
p  Ip 

non-incident  branch  /  incident  branch 


D3050-14 

Left  (R/L)  Right  (L/R) 


eb 

0.3842 

e2 

0.6118 

ec 

0.4451 

e3 

0.6278 

e4 

0.6185 

Average 

0.41465 

0.6194 

STDEV 

0.043063 

0.0080 

105 


Left  Side  Launched 


a) 

Right  Side  Launched 


b) 

Figure  6.46  Measured  data  and  corrected  data  taken  from  the  directional  couplers  of  a 
single  sample.  The  values  are  in  the  P2/P1  ratio  form. 
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<  2.79  jLim  ► 


* -  3.53  pm  * 

Figure  6.47  Detailed  dimensions  of  a  directional  coupler  sample 

In  order  to  compare  the  data  with  the  theoretical  results,  precise  measurements  of  the 
waveguides  are  needed.  As  it  can  be  seen  in  Figure  6.15  the  waveguides  constructed  are 
not  perfectly  square.  Also,  the  waveguide  width  is  not  identical  to  the  template  mask. 
Measurements  were  taken  to  determine  the  exact  shape  of  the  etched  waveguide.  These 
measurements  indicated  that  instead  of  the  waveguide  being  a  square  it  is  actually 
trapezoidal  in  shape.  Figure  6.47  shows  the  results  of  the  measurements,  the  grayed  area 
is  the  GaN  core  region  and  the  white  area  is  the  AlGaN  cladding  region.  With  these 


Table  6.  Widths  of  the  directional  coupler.  Height  is  2.09  pm.  Sample  number  is  A 12941 


Width  (nm) 


Coupler  # 

coupling 

region 

Single 

waveguide 

#4 

2.79/3.53 

1.48/2.14 

#5 

2.89/3.6 

#17 

3.03/3.74 

#28 

3.08/3.8 
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Eperimental  Data  with  theoretical  prediction  for  A12941 


Coupling  length  (pm) 


Figure  6.48  Comparison  of  the  data  from  sample  A12941  with  theoretical  predictions  for 
1.71  pm  width  branches. 

results  one  can  calculate  the  slope  angle  0  to  be  79.72°  (  ~  80° ),  the  upper  width  of  the 
core,  xi,  is  3.15  pm,  and  the  lower  width  of  the  core,  X2,  is  3.44  pm.  Table  6  shows  the 
measurements  of  other  directional  couplers  in  the  same  sample 
With  these  measurements,  the  averages  were  then  plotted  and  compared  with  the 
theoretical  results  from  Marcatili’s  method  and  the  beam  propagation  method.  Since  the 
waveguide  is  trapezoidal  an  effective  width  can  be  found  that  would  match  the  behavior 
of  the  directional  coupler.  Figure  6.48  shows  the  corrected  experimental  data  and  the 
plotted  theoretical  behavior.  For  Marcatili’s  method,  the  sum  of  square  error  was  done  for 
various  coupling  lengths.  The  result  was  a  width  of  ~3.42  pm  for  the  coupling  region, 
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corresponding  to  a  coupling  length  of  -290  pm,  or  -1.71  pm  for  the  single  waveguide 
width.  This  width  is  within  the  boundaries  of  the  core  width  calculated  from  Figure  6.47. 
With  the  data  point  labeled  2  the  SSE  (sum  of  square  errors)  between  the  Marcatili’s 
results  and  the  data  is  0.222  (  with  a  coupling  length  of  289.4  pm  ).  Without  this  data 
point  the  SSE  was  reduced  by  half  to  0.1 1 1  (  and  the  coupling  length  was  290.021  pm  ). 
The  structures  with  a  normalized  intensity  above  0.5  showed  a  crossover  in  intensity, 
when  the  maximum  intensity  switches  from  one  branch  to  the  other.  The  data  didn’t 
support  the  beam  propagation  method  simulation,  which  predicted  a  much  smaller 
coupling  length  for  the  same  1.71pm  widths  of  the  branches,  as  illustrated  in  Figure  6.48. 
The  width  in  the  beam  propagation  simulation  would  have  to  be  much  larger  than  the 
actual  physical  device  to  obtain  the  same  coupling  length.  The  SSE  for  the  beam 
propagation  method  was  0.7  without  data  point  2. 

A  second  sample  was  analyzed  in  the  same  fashion  as  sample  A 12941.  This  time  a 
different  set  of  coupling  lengths  were  used.  Table  7  shows  the  dimensions  on  some  of  the 
structures  in  the  sample.  Figure  6.49  shows  the  corrected  data  with  the  Marcatili  and  the 


Table  7.  Widths  of  the  directional  coupler.  Height  is  1.986  pm.  Sample  number  is 
A28542 


Width  (nm) 


Coupler  # 
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3.12/3.83 
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#21 
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#23 

3.17/3.87 
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Experimental  data  and  theoretical  results  for  a28542 


Coupling  length  (pm) 


Figure  6.49  Comparison  of  the  data  from  sample  A28542  with  theoretical  predictions  for 
1.73  pm  width  branches 

beam  propagation  methods.  This  particular  sample  had  more  experimental  errors  but  it 
still  tracked  the  Marcatili  simulation  better  than  the  beam  propagation  simulation.  The 
SSE  for  the  Marcatili  simulation  was  0.154  and  for  the  beam  propagation  simulation  is 
4.135. 

With  this  mask  set,  the  processing  for  the  active  directional  coupler  was  also  developed. 
Depositing  the  contact  on  top  of  the  waveguide  was  difficult  but  achievable.  The  “Top 
Contact”  process  in  APPENDIX  C.2  shows  the  process  followed  to  deposit  the  ohmic 
contact  on  the  directional  coupler.  The  composition  of  the  contact  is  Ti/Al/Ti/Au  and  the 
heights  used  for  each  are  detailed  in  the  process  follower.  Figure  6.50  shows  an  SEM 
image  of  the  contact  on  the  directional  coupler.  Notice  that  although  the  contact  strip  is 
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Figure  6.50  SEM  image  of  contact  on  top  of  the  directional  coupler. 


up  to  the  edge  it  didn’t  flow  to  the  sidewall.  The  trench  and  the  metal  contacts  processes, 
see  Figure  6.34,  were  created  with  little  difficulty,  see  APPENDIX  C.2.  In  order  to 
activate  the  ohmic  contact,  annealing  is  required. 

A  first  annealing  attempt  is  shown  in  Figure  6.51.  The  annealing  was  done  in  a  rapid 
thermal  annealer  (  RTA  )with  a  nitrogen  flow.  Figure  6.51  a)  and  b)  were  obtained  from 
the  same  directional  coupler  at  different  locations.  Note  that  Figure  6.51  a)  seems  to  have 
a  different  morphology  than  b).  Figure  6.52  seems  to  indicate  that  the  difference  in 
morphology  is  due  to  bottom  AlGaN  contact  acting  as  a  radiator  and  subjecting  the  areas 


a)  b) 

Figure  6.51.  First  annealing  attempt  with  an  RTA  at  850°  C  for  30s 


Ill 


Figure  6.52  Lower  magnification  than  figure  Figure  6.51 
surrounding  it  to  less  thermal  stress.  Additional  examination  under  the  SEM  showed  little 
surface  difference.  The  SEM  images  of  Figure  6.50  were  taken  from  these  two  regions. 
Figure  6.53  is  inconclusive  as  to  the  cause  of  the  difference  in  appearance.  The  difference 
could  have  been  from  oxygen,  nitrogen,  or  a  secondary  Ti  peak. 

The  first  attempt  at  creating  a  bridge  was  not  very  successful.  The  composition  used 
was  Ti/Au  with  the  thickness  of  200/500  A.  Altering  the  bridge  composition  to 


0.80  1.60  2.40  3.20  4.00  4.80  5.60 

Figure  6.53  EDAX  analysis  of  the  contacts  after  annealing.  The  Black  outline  is  from  the 
center  section  and  the  solid  graph  is  from  the  morphologically  different  tip. 
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Figure  6.54  SEM  images  of  the  bridges 


300/4000A  permitted  us  to  obtain  nearly  100  percent  survivability.  In  APPENDIX  C.2 
starting  with  the  process  “Contact  Opening  and  Reflow”  and  ending  with  the  process 
“Bridge”  give  a  description  of  how  the  bridges  were  created  and  their  compositions. 
Figure  6.54  shows  the  SEM  images  of  the  bridges. 

The  I-V  curve  of  the  device  is  shown  on  Figure  6.55.  The  device  resistance  was 
~  12.3  Os.  A  low  resistance  was  expected,  since  in-house  efforts  were  not  producing 


l-V  curve  of  Sample  A235-4-1 


Figure  6.55.  I-V  curve  of  active  device 
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semi-insulating  GaN  layer  at  the  time.  When  ~5  V  was  applied  to  an  active  device  it 
showed  an  unusual  behavior.  The  optical  signal  was  coming  out  of  both  branches.  When 
a  voltage  was  applied  both,  turned  off.  One  of  the  samples  cracked  right  through  the 
waveguide,  see  Figure  6.56.  The  circled  region  in  Figure  6.56  is  where  the  electrical 
probe  was  making  contact.  The  sample  broke  right  through  the  location  where  the  other 
probe  was  located. 

Another  sample  was  tested  and  the  probes  were  observed  to  be  moving,  thermal 
expansion  was  suspected.  When  the  sample  was  taken  out  of  the  holder  the  explanation 
surfaced.  The  sample  was  getting  so  hot  that  it  melted  a  plastic  holder.  A  new  holder  was 
made  out  of  aluminum  and  installed.  The  sample  holder  was  insulated  from  the  rest  of  the 
set  up.  When  voltage  (<=  6V)  was  applied,  small  “sparks”  were  seen  to  be  coming  out  of 


Figure  6.56.  Active  directional  coupler  that  cracked  after  applying  voltage. 


114 


the  negative  probe  and  damaged  the  contact  pad  underneath.  The  probes  didn’t  “move” 
this  time.  This  showed  that  the  plastic  holder  being  melted  was  the  cause  of  the  probe 
movement.  An  approach  that  was  later  attempted  was  to  insulate  the  contact  pads  with 
SiN.  This  ensured  the  functionality  of  the  device  as  a  field  only  device. 

During  the  course  of  the  manufacturing  of  the  device  structures,  slight  variations  were 
noted  on  the  devices  from  one  side  to  the  other.  Figure  6.57  shows  a  view  in  the  mask 
template  near  the  junction  of  the  directional  coupler.  There  is  a  slight  variation  in  the 
width  of  the  two  branches.  This  error  in  the  mask  could  have  been  a  mayor  contributor  of 


Figure  6.57  Junction  in  a  directional  coupler  directly  on  the  mask 
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Figure  6.58  Transferred  pattern  into  the  AlGaN/GaN/AlGaN  material 


the  uneven  symmetry  of  the  data  seen  in  Figure  6.46.  Figure  6.58  shows  the  pattern 
transferred  to  the  material,  after  ICP  etching.  The  difference  between  the  branches  is 
nearly  unnoticeable  at  this  stage.  Although,  they  could  have  been  enough  for  the 
discrepancy  in  the  data. 

Third  Mask  Set  Processing,  and  Results 

Another  mask  set  was  needed  in  which  the  errors  in  this  mask  set  was  corrected.  With 
the  experience  gained,  modifications  were  made  to  the  mask  to  make  the  alignment 
easier.  The  mask  arrangement  is  given  in  APPENDIX  C.2.  When  the  new  mask  set  was 
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sent  to  the  company,  and  after  informing  them  of  the  errors  that  was  encountered  they 
were  able  to  fix  the  error  in  their  machine.  This  time,  along  with  fixing  software  errors, 
the  new  mask  set  was  done  with  finer  resolution.  The  new  resolution  was  0.2  pm.  This 
brought  back  the  problems  in  the  junction,  see  Figure  6.36  and  Figure  6.37.  This  was 
because  the  thin  feature,  circled  in  these  figures,  were  smaller  than  in  the  previous  mask 
set  template. 

In  order  to  resolved  it  this  time,  a  thinner  photoresist  was  used.  The  new  process  is 
listed  in  APPENDIX  C.3.  These  process  followers  were  arrived  at  after  several  iterations. 
With  the  new  mask  set  and  process  followers,  a  much  finer  detail  was  achieved.  Figure 
6.59  shows  the  junction  on  the  directional  coupler.  This  time,  with  the  new  mask  set,  both 


Figure  6.59  Junction  of  a  directional  coupler  with  the  new  mask  template  and  the  process 
outlined  in  APPENDIX  C.3 
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Figure  6.60  Different  junction  done  with  the  third  mask  template 

branches  are  identical  in  width.  In  Figure  6.60  the  same  mask  set  was  used  to  create  the 
pattern.  A  discrepancy  is  seen  between  the  branches,  Figure  6.60  inside  the  brackets.  This 
is  also  present  in  the  mask  template  and  is  due  to  quantization  error  when  transferring  the 
pattern  from  the  software  into  the  mask  template. 

Since  the  AlGaN  couldn’t  be  made  insulating  enough  for  the  device  to  work  properly, 
see  previous  section,  a  work  around  was  devised.  This  time  the  contacts  were  insulated 
from  the  AlGaN/GaN  material  by  a  SiN  layer,  as  illustrated  in  Figure  6.62.  In  that 
configuration,  the  electric  field  that  the  core  would  be  subjected  to  would  be  a  lot  smaller. 
This  is  because  the  voltage  is  not  concentrated  in  the  core  but  distributed  between  the 
contacts.  Therefore,  a  larger  voltage  is  required  to  operate  the  directional  coupler. 
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Figure  6.62  The  targeted  new  configuration  for  the  active  device 


A  new  sample  was  prepared  with  this  configuration  using  the  new  mask  template  and 
the  process  detailed  in  APPENDIX  C.3.  A  uniform  layer  of  1000  A  of  Si3N4  was 
deposited  on  the  sample,  with  a  new  PECVD  system.  After  the  metal  contacts  were 
deposited,  the  sample  was  placed  on  an  RIE  chamber  to  etch  the  Si3N4  where  it  was 
exposed.  This  should  have  left  only  the  Si3N4  under  the  contacts.  The  bridge  was  then 


Comparison  of  data  collected  with  corrected  data  for  A1 8711 
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Figure  6.61  Measured  data  and  corrected  data  taken  from  the  directional  couplers  of  a 
sample.  The  values  are  in  the  P2/P1  ratio  form. 
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made  and  the  ends  polished.  Data  was  taken  on  the  devices  before  applying  a  voltage. 
Figure  6.61  shows  the  data  gathered  from  the  devices.  As  it  can  be  seen  the  measured  and 
the  corrected  data  vary  very  little  from  each  other,  in  contrast  with  that  of  Figure  6.46. 
This  can  be  attributed  to  the  improved  mask  template,  which  has  finer  resolution  and  has 
corrected  the  previous  errors. 

In  examining  the  data,  discrepancies  started  to  show.  The  width  of  the  coupling  section 
was  measured  with  in  the  SEM  and  turned  out  to  be  3.55/4.21  pms  with  a  height  of 
2.42  pm,  see  Figure  6.47.  This  would  give  us  a  single  waveguide  width  between  1.91- 
2.02  pm  The  data  pointed  to  a  smaller  coupling  length  than  Marcatili’s  approach  would 
indicate  for  these  dimensions,  see  Figure  6.48  and  Figure  6.49.  Using  Equations  (  4-10  ) 


Experimental  data  with  LSE  fit  of  power  equations  for  A2871 1 


Coupling  length  (pm) 


Figure  6.63  Least  square  fit  of  the  data  with  Equations  (  4-10  )  and  (4-11  ).  Coupling 
length  is  63.954  pm 
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Comparison  of  data  and  BP  with  air  cover  for  a2871 1 


Coupling  length  (pm) 


Figure  6.64  Comparison  of  the  data  from  sample  A287-1-1  and  the  beam  propagation 
method  with  air  covering  the  sample. 

and  (4-11  )  with  a  =  0  and  doing  a  least  square  fit,  one  arrives  at  a  coupling  length  of  ~ 
64  pm,  the  sum  of  square  errors  is  0.014  for  this  value,  see  Figure  6.63.  This  would  be  the 
approximate  coupling  length  that  we  would  have  arrived  for  a  1.71-1.73  pm  branches, 
with  a  3.42  -  3.46  pm  total  width  of  the  coupling  region,  using  the  beam  propagation 
method,  see  Figure  6.64.  Investigation  of  the  surface  of  the  sample  pointed  to  the  reason 
of  the  discrepancy  in  the  data. 

Examining  the  surface  of  the  sample  using  the  EDAX  tool  pointed  to  the  cause  of  the 
discrepancy.  Figure  6.65  illustrates  that  Si3N4  is  still  on  the  surface.  The  index  of 
refraction  of  the  grown  Si3N4  is  ~  2.26,  which  is  very  close  to  that  of  the  core  ~  2.304. 
This  proximity  to  the  index  of  refraction  of  the  core  will  affect  the  coupling  length  of  the 
directional  coupler.  Attempts  were  made  to  determine  the  height  of  the  Si3N4  left  on  the 
sample.  Figure  6.66  shows  the  cross  sectional  view  of  a  pad  after  milling.  From  the  cross 
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Figure  6.65  ED  AX  of  a  bare  spot  on  the  A2871 1  sample 
sectional  view,  the  Si3N4  could  not  be  discerned,  making  the  height  impossible  to  be 
determined.  An  EDAX  done  on  top  of  the  metal  contact  showed  that  Si3N4  was  present 
along  with  the  metal. 

A  series  of  simulations  were  done  on  the  beam  propagation  method  and  the  Marcatili 
method  in  which  the  entire  covering  material  was  changed  from  air  to  Si3N4.  The 
Marcatili  method  gave  a  coupling  length  of  ~  87  pm  for  branches  width  of  1.7  pm,  which 


Figure  6.66  Cross-sectional  view  at  a  contact  pad  for  sample  A287-1-1 
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Comparison  of  data  and  BP  model  for  a28711 


Coupling  length  (pm) 


Figure  6.67.  Comparison  of  the  data  and  simulations  with  the  beam  propagation  method. 

is  bigger  than  the  suggested  by  the  data,  see  Figure  6.63.  With  the  beam  propagation 
simulation  two  approaches  were  attempted  one  in  which  Si3N4  covered  the  entire 
structure  and  one  in  which  only  a  0.1  pm  fdm  covered  the  structure.  This  thickness  was 
used  because  it  was  the  maximum  thickness  that  the  sample  could  have,  assuming  that  the 
RIE  with  freon  had  little  or  no  effect  on  the  Si3N4.  The  simulations,  using  the  beam 
propagation  method,  seems  to  have  additional  modes  present  as  seen  in  Figure  6.67. 

Voltage  was  applied  to  the  structure  to  attempt  to  see  if  the  optical  signals  would  switch 
or  modulate.  After  the  voltage  was  applied  there  was  no  apparent  change  in  the  intensity 
of  the  output  signal.  An  I-V  curve  of  one  of  the  structure  was  obtained  and  a  possible 
reason  emerged.  Figure  6.68  shows  the  I-V  curve  of  the  structure.  The  Si3N4  was  not 
acting  as  an  insulator.  When  the  sweeping  voltage  was  set  to  ±10  V  the  structure  was 
acting  as  a  diode  with  a  forward  a  turn  on  voltage  of  ~  0.7  V.  When  the  sweep  was 
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l-V  curves  of  structure  in  sample  A28711 


Figure  6.68  I-V  curve  of  a  structure  in  sample  A287-1-1.  The  on  resistance  is  between 
10-12  k^s 

increased  beyond  15  V  the  device  appeared  to  have  gone  beyond  breakdown  and  from 
then  on  it  went  into  a  resistive  behavior,  no  matter  how  small  the  voltage  sweep  was  set 
at.  While  taking  the  I-V  curves,  and  after  device  breakdown,  the  curves  were  jittery,  the 
slope  and  shape  of  the  curve  was  oscillating.  It  was  suggested  that  there  might  be  some 
type  of  traps  in  the  structure  accounting  for  the  jittery  behavior.  Apparently  the  structure’s 
resistive  behavior  prevented  the  electric  field  in  the  core  region  to  be  sufficient  to 
noticeably  change  the  refraction  index. 

In-house  efforts  to  obtain  semi-insulating  (  compensated  )  and  p-type  doping  have  been 
improving.  A  sample  was  prepared  in  an  attempt  to  make  the  GaN  layer  insulating,  see 
Figure  6.6.  In  this  structure  the  cladding,  Al0.05Ga0.95N,  layers  were  both  doped  n-type,  ~ 
1018  cm'2.  For  this  purpose  Mg  was  incorporated  in  the  growth  to  compensate  the 
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l-V  curve  for  the  structure  in  sample  A421-3-1 


Figure  6.69  I-V  curve  of  sample  A421-3-1  in  which  attempts  were  made  to  obtain 
insulating  GaN 

unintentionally  doped  n  concentration.  Si3N4  was  deposited  prior  to  the  contacts  on  the 
entire  surface,  ~  1000  A.  This  time  no  attempts  to  remove  it  were  made.  For  the  contacts, 
Ti/Au  was  used,  since  this  has  been  shown  to  be  resistive  for  n-type  AlGaN  without 
alloyng.37 

Testing  was  started  on  the  sample.  Prior  to  testing  the  structure  in  the  optical  set  up  an 
I-V  curve  was  made  of  the  structure.  It  was  a  different  device  than  the  one  to  be  tested. 
Figure  6.69  shows  the  I-V  curves  taken  from  a  sample.  From  this  figure  it  is  seen  that, 
after  breakdown,  the  device  goes  into  a  resistive  mode  with  a  resistance  of  ~5  kfl  The 
current  of  the  device  was  monitored  during  testing,  in  the  optical  set  up,  to  prevent  any 
degradation  of  the  structure,  see  Figure  6.17.  The  structure  presented  some  diode 
behavior  during  testing,  with  one  polarity  of  the  voltage  presenting  slightly  less  resistance 
than  the  other.  The  structure  should  have  presented  a  highly  resistive  behavior. 


a)  Lights  off  b)  Light  on 

Figure  6.70  Image  obtained  from  the  output  of  a  directional  coupler  of  sample  A421-3-1. 

Data  was  collected  on  one  of  the  structures  in  the  sample.  Figure  6.70  shows  the  output 
of  the  directional  coupler.  Behavior  as  a  directional  coupler  was  immediately  noticed. 
Light  was  launched  in  the  left  branch,  the  opposite  branch  from  that  shown  in  Figure 
6.70,  and  light  came  out  in  the  left  branch.  Once  all  the  data  was  collected  from  the 
device,  voltage  was  applied  and  increased  slowly.  During  that  time,  the  current  was 
monitored  for  a  possible  breakdown  of  the  device.  The  voltage  increase  was  stopped 
when  the  current  started  to  increase  rapidly.  The  values  of  the  current  and  voltage  were 
0.76  pA  and  25  V  respectively.  The  intensity  of  the  output  appeared  to  change.  The 
voltage  was  turned  off  and  on  to  verify  a  change  in  intensity.  Figure  6.71  shows  instances 
in  which  the  voltage  was  turned  on  and  off.  The  difference  in  intensity  is  more  noticeable 
in  the  intensity  profile,  section  c)  of  the  figure.  There  was  no  apparent  change  in  the 
branch  that  had  no  optical  signal  in  the  output.  A  possible  reason  is  that  the  branch  is 
below  the  noise  level  and  any  variation  would  be  difficult  or  impossible  to  detect. 
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a)  Voltage  off  b)  Voltage  on 


1  21  41  61  81  101  121  141  161  181  201  221  241  261  281 

Sample  # 

c)  Intensity  profile 

Figure  6.71  Light  intensity  variation  with  the  application  of  25  V  on  a  device  in  sample 
A421-3-1. 

The  output  of  the  camera  was  manipulated  to  optimize  the  contrast  and  intensity.  That 
was  done  electronically  rather  than  with  image  manipulation.  This  accentuated  the 
variation  of  the  signal,  by  removing  a  portion  of  the  dc  component  of  the  image  and 
focusing  on  the  intensity  variation.  Figure  6.72  shows  the  results  of  the  modulation  of  the 
directional  coupler.  The  noise  in  the  signal  is  due  to  building  vibrations. 


Intensity  (  a.  u.  and  dc  offset) 
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Directional  coupler  modulation 


Figure  6.72  Modulation  of  the  directional  coupler.  Image  enhancements  were  done  to 
accentuate  the  intensity  variation. 


CHAPTER  VII 


DISCUSSION  OF  RESULTS  AND  CONCLUSION 


Developing  the  process  for  etching  the  AlGaN/GaN  material  system  proved  to  be  more 
difficult  than  initially  anticipated.  This  was  due  mostly  to  the  newness  of  the  material  and 
the  limited  experience  in  processing  it.  This  is  illustrated  by  the  initial  search  for  an 
appropriate  masking  material  and  the  process  development  as  seen  in  appendices  B  and 
C.  Since  GaN/AlGaN  is  highly  tolerant  of  harsh  environments,  the  normal  masking 
materials  that  would  fare  well  on  other  materials  would  only  be  marginal  at  best.  The 
photoresist,  normally  used  as  a  masking  material,  etches  a  little  faster  than  the  AlGaN 
material  system.  This  would  be  adequate  for  shallow  structures,  but  not  very  good  for 
deep  structures,  such  as  a  directional  coupler.  Figure  7.1  shows  photoresist  after  being 


a)  b) 


Figure  7.1  Photoresist  before  metal  mask  deposit,  a)  Junction  region  of  a  directional 
coupler,  b)  comer  of  a  letter 
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processed  into  a  mask.  The  sidewalls  of  the  photoresist  shown  in  Figure  7.1b  are  not  very 
promising  for  a  mask  material.  A  problem  that  would  immediately  surface  would  be  the 
sidewall  slope  and  roughness.  For  the  targeted  structure  of  this  work  the  sidewalls  needed 
to  be  steep  and  as  smooth  as  possible.  Finding  Ni  as  a  mask  was  very  advantageous.  Only 
1000  A  of  material  was  needed  to  etch  ~  2  pm  of  material,  it  also  permitted  steep 
sidewalls.  Placing  the  electrode  on  top  of  the  directional  coupler  presented  an  alignment 
challenge.  The  electrodes  were  less  than  3  pm  wide  and  ~3  mm  long  and  they  required  a 
laborious  alignment  procedure.  Figure  7.2  shows  the  successful  placement  of  the  contact 
on  top  of  the  coupling  section  of  a  directional  coupler. 

The  structure  that  this  research  focused  on  is  shown  in  Figure  7.3.  This  is  a  complicated 
structure  to  simulate.  As  seen  in  Figure  7.4,  in  the  vertical  direction  the  evanescent  waves 


Figure  7.2.  Metal  contact  on  top  of  the  coupling  region  of  a  directional  coupler. 
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Figure  7.3  Side  view  of  the  single  mode  waveguide  structure  for  AlGaN 

extend  into  the  AlGaN  cladding  region.  In  a  two-dimension  structure  this  would  be  close 
to  the  cut  off  region.  In  the  horizontal  direction  the  evanescent  waves  don’t  penetrate  the 
cladding  region,  which  is  air  in  that  direction.  That  would  be  far  from  cut-off  in  a  two 
dimensional  structure.  Having  both  of  these  conditions  present  is  what  makes  this 
structure  complicated  to  simulate. 

The  Marcatili  method  is  well  suited  for  simulation  in  the  horizontal  direction,  where  air 
(  n  =  1)  acts  as  the  cladding  region.  In  that  direction,  the  index  of  refraction  has  abrupt 


X  (|xm) 


Figure  7.4  Intensity  profile  in  the  single  mode  waveguide 
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changes  and  the  Marcatili  method  is  tailored  for  abrupt  changes,  see  CHAPTER  2.  But 
for  the  parts  that  correspond  to  the  shaded  area  in  Figure  2.4,  where  AlGaN  is  present,  its 
suitability  could  be  sketchy.  The  beam  propagation  method  is  more  suited  for  the  vertical 
direction  because  the  change  in  the  index  of  refraction  is  smaller  and  the  BPM  is  better 
with  small  refraction  index  increments,  see  CHAPTER  2. 

In  the  1.4  pm  waveguide  width  simulations,  Figure  6.25  and  Figure  6.26,  there  was  a 
discrepancy  between  the  methods.  The  Marcatili  method  showed  a  waveguide  that 
shouldn’t  be  guiding,  Figure  6.26,  while  the  beam  propagation  method  showed  a  guiding 
waveguide.  Although  the  mask  template  had  this  dimension,  the  actual  waveguide  was 
wider  than  this.  The  thinnest  waveguide  tested  was  ~  1.71  pm.  In  Figure  7.5  both  the 
BPM  and  the  Marcatili  method  are  compared  side  by  side.  While  the  Marcatili  model  is 
very  close  to  cut  off,  the  BPM  is  guiding  very  well.  Figure  7.6  shows  the  output  of  the 
1.71  pm  branch  width  directional  coupler.  It  points  toward  the  beam  propagation  method 
as  a  better  representation  of  the  wave  being  guided  in  the  waveguide.  This  can  be 
explained  by  knowing  that  one  of  the  assumptions  of  the  Marcatili  method  is  that  the 


Power  in  Waveguide  Mode 


Propagation  Direction  (pm) 


B-V  Curve 


a)  Beam  propagation  b)  Marcatili 

Figure  7.5  Comparison  of  the  1.71pm  waveguides. 


Figure  7.6.  Output  from  the  1.71  pm  branch  width  directional  coupler,  sample  A 12941 


optical  signal  needs  to  be  far  from  cut  off,  see  CHAPTER  2.  When  the  optical  signal  is 
close  to  cut  off  in  the  waveguide,  the  Marcatili  method  will  not  accurately  account  for  the 
behavior  of  the  optical  signal  in  the  waveguide,  see  reference  22.  Figure  7.5b  shows  a 
waveguide  too  close  to  cut  off  to  be  effectively  propagating  light.  The  beam  propagation 
method  shows  an  optical  signal  that  is  propagating  effectively  through  the  waveguide, 
giving  a  better  representation  of  what  it  is  observed.  Also,  the  errors  incurred  in  the 
paraxial  approximation  were  not  a  major  factor  in  the  behavior  with  a  single  waveguide 
simulation. 

For  the  directional  coupler  the  results  were  different.  Figure  6.48  and  Figure  6.49  favors 
the  Marcatili  approximation  for  establishing  the  coupling  length  in  the  directional 
coupler.  The  paraxial  approximation  of  the  beam  propagation  method  neglects  large 
angles  from  the  z-direction  and  thus  introduces  an  error  in  the  calculation  of  the  coupling 
length.  A  directional  coupler  is  very  sensitive  to  these  errors  and  thus  gives  an  erroneous 
representation  of  directional  coupler  behavior.  In  a  different  way,  the  beam  propagation 


133 


B-V  Curve 


Figure  7.7  Mode  chart  of  the  interaction  region  for  a  1.71  gin  branch  width  directional 
coupler 

method  assumes  small  variations  in  the  index  of  refraction  not  only  in  the  propagation 
direction,  but  also  in  the  transverse  direction.  The  abrupt  changes  in  the  air-core  interface 
introduce  errors  in  the  propagation  constants  of  the  two  modes,  odd  and  even,  in  the 
coupling  region.  This  in  turn  gives  us  an  erroneous  coupling  length.  The  B-V  curves  from 
the  Marcatili  method,  shown  in  Figure  7.7,  shows  one  of  the  modes  near  cut  off.  In  spite 
of  the  method  of  being  unreliable  near  mode  cut-off  its  results  were  close  to  what  was 
expected,  see  Figure  6.47  and  surrounding  text.  The  reasoning  for  this  is  that  the  coupling 
length  is  determined  by  the  difference  between  the  modes,  and  the  errors  in  the  Marcatili 
approach  had  little  effect  on  the  difference  for  this  case. 

When  the  medium  surrounding  the  directional  coupler  has  a  smaller  step  increase  in  the 
index  of  refraction,  the  Marcatili  method  will  no  longer  able  to  predict  the  behavior  of  the 
waveguide.  That  is  because  the  shaded  areas,  see  Figure  2.4  and  Figure  6.3,  can  no  longer 
be  neglected.  When  SiN  was  introduced  surrounding  the  structure,  it  was  enough  to  be  in 
a  shaded  region  where  neither  the  beam  propagation  nor  the  Marcatili  method  could 
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Figure  7.8.  Structure  of  the  modulating  directional  coupler 


establish  theoretically  the  coupling  length.  Although  this  time,  the  beam  propagation 
method  was  much  closer  than  before.  Possible  variations  in  the  Al%,  the  core  size,  and 
the  index  of  refraction  could  have  also  accounted  for  this  discrepancy.  Though  it  is  still 
believed  that  better  modeling  methods  are  needed  in  this  region. 

Modulation  was  obtained  in  the  structure  of  Figure  7.8.  That  indicated  that,  indeed  the 
electo-optic  effect  is  present  in  the  material,  although,  confirmation  of  its  value  was  not 
possible.  The  modulation  depth  was  much  lower  than  was  originally  sought  after.  In  this 
structure,  after  device  breakdown,  SiN  didn’t  act  as  an  insulator,  as  expected.  Another 
unexpected  result  was  the  entire  structure  having  a  low  resistive  value  ~5  kfi  after  SiN 
breakdown,  see  Figure  6.69.  In  other  reaserch,  this  dielectric  has  been  shown  to  fail  on 
III-V  materials.37  Because  Mg  diffuses,  it  was  added  to  the  growth  for  only  0.5  pm  of  the 
GaN  region,  see  Figure  7.8,  the  rest  of  the  GaN  was  undoped.  This  was  done  in  the 
expectation  that  Mg  would  diffuse  over  the  entire  GaN  region  and  create  a  compensated 
insulator.  The  weak  modulation  could  be  attributed  to  the  fact  that  the  voltage  is  applied 
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to  a  distance  of  more  than  50  times  as  great  as  intended,  >  50  pm  instead  of  0.8  pm.  That 
reduced  the  electric  field  intensity  in  the  core,  thus  we  have  only  a  small  index  variation 

Future  Work 

Through  out  this  project  one  thing  was  constantly  surfacing,  the  material  is  still  in  the 
development  stage.  The  material  growth  needs  to  improve  the  development  of  insulating, 
or  semi-insulating,  layers  with  out  requiring  compensation.  Better  control  and  improved 
p-type  layers  are  needed.  Additionally,  a  reduction  in  defect  density  should  be  pursued 
since  this  could  reduce  the  background  n  concentration.  In  the  modeling  for  the  structure 
under  investigation,  more  robust  models  should  be  pursued.  Alternatives  include  the  use 
of  the  finite  element  approach  to  obtain  better  B-V  curves  than  can  be  obtained  by  the 
Marcatili  method.  The  beam  propagation  method  is  a  good  tool  to  use  to  model  the 
behavior  of  the  entire  structure.  This  includes  the  study  of  small  imperfections  that  can  be 
added  to  the  structure,  which  could  be  produced  by  the  manufacturing  process,  and 
observe  the  guides  behavior.  Unfortunately,  its  domain  of  operation  does  not  include  the 
abrupt  changes  in  the  index  of  refraction  in  the  transverse  direction.  If  possible, 
modifications  should  be  sought  to  make  it  less  susceptible  to  abrupt  index  changes 
without  sacrificing  much  of  its  computational  advantages.  Another  possibility  is  to  use 
the  finite  difference-time  domain  (  FD-TD  )  approach  for  the  full  structural  analysis. 
The  problem  is  that  this  approach  is  very  computational  intensive  and  memory  hungry. 
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CONFIGURATION  OF  MASK  1 
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APPENDIX  B 


ETCHING  MASK  AND  ICP  RECIPE  DEVELOPMENT 


Recipe  Number 

ICP  Recipe  detail 

1 

(Paul’s  Recipe) 

Cl2 

4  seem 

BC13 

15  seem 

Pressure 

5  mTorr 

RF  Power 

75  W 

ICP  Power 

300  W 

2 

(Dr  Shull’s  Recipe) 

Cl2 

32  seem 

BC13 

8  seem 

Ar 

5  seem 

Pressure 

2  mTorr 

RF  Power 

150  W 

ICP  Power 

500  W 

3 

Cl2 

20  seem 

Ar 

5  seem 

Pressure 

5  mTorr 

RF  Power 

150  W 

ICP  Power 

500  W 

4 

Cl2 

20  seem 

BC13 

20  seem 

Ar 

5  seem 

Pressure 

2  mTorr 

RF  Power 

150  W 

ICP  Power 

500  W 
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On  the  sample  there  were  cluttered  (dirty) 
areas  and  clean  areas. 
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B.2  ICP  Recipe  Development  for  Sidewall  Improvement 
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APPENDIX  C 


PROCESS  FOLLOWERS 
C.l  First  process  follower  for  passive  Waveguides 


Step 

Process  Step 

1 

Pre  SiN  SAMPLE  PREPARATION  (OPTIONAL): 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Prepare  surface  with  a  dip  of  Ammonia  Hydroxide:DI  (1:20)  for  1  min 

□  Deposit  1000  A  of  SiN 

2 

SAMPLE  PREPARATION: 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

3 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

4 

PMGI  COAT: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

5 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

6 

EXPOSE  1813  with  BASE  MASK: 

(Note  which  exposure  was  done  first) 

□  Align 

□  4.0  sec  Exposure  (405  nm)  on  Mask  Aligner  @  20mW/cm2  optical  intensity 

7 

1813  DEVELOP: 

□  40  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

8 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

9 

DEEP  UV  (  254  nm  ): 

□  100  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 
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Step 

Process  Step 

10 

PMGI  DEVELOP: 

□  50  sec  101  Develop  Immersion 

□  4X  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

11 

DEEP  UV  (  254  nm  ): 

□  100  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

12 

PMGI  DEVELOP: 

□  50  sec  101  Develop  Immersion 

□  4X  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

13 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

14 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

15 

EVAPORATE  BASE  METAL: 

□  Evaporate  1000  A  Ni 

16 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Remove  base  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

17 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

18 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

19 

SiN  DEFINITION  (OPTIONAL): 

□  Etch  using  the  dual  chamber  RIE  reactor 

20 

ETCH  AlGaN: 

□  Etch  using  the  ICP  reactor 

21 

PREPARATIONS: 

□  Mix  (1:1)  BOE:DI  H2O  using  at  least  25  ml  BOE  with  25  ml  DI  H2O  in  proper 
container 

22 

REMOVE  SiN  (OPTIONAL): 

□  3  min  BOE:DI(l:l) 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

C.2  Process  follower  for  second  Mask  set 


Purpose  Mask  Process 

Active  waveguide  Nitride  etch  (dark  field)  WG  definition 


Step. 

Process  Step 

1 

Pre  SiN  SAMPLE  PREPARATION  (OPTIONAL): 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  1  min  ammonium  hydroxide 

□  Prepare  surface  with  Ammonia  Hydroxide:DI  (1:20)  for  20  sec 

□  Deposit  1000  A  of  SiN 

2 

SAMPLE  PREPARATION: 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

3 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

4 

PMGI  COAT: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

5 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

6 

EXPOSE  1813  with  Contact  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (8.60  settings) 

@  20mW/cm2  optical  intensity 

7 

1813  DEVELOP: 

□  60  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 
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Step. 

Process  Step 

8 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

9 

DEEP  UV  (  254  nm  ): 

□  100  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

10 

PMGI  DEVELOP: 

□  60  sec  SAL  101  Develop  Immersion 

□  4x  bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

11 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

12 

DEEP  UV  (  254  nm  ): 

□  100  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

13 

PMGI  DEVELOP: 

□  60  sec  SAL  101  Develop  Immersion 

□  4x  bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

14 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

15 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

16 

EVAPORATE  Contact  METAL: 

□  Evaporate  1000  A  Ni 

17 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

18 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 
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Step. _ Process  Step _ 

1 9  02  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

20  SiN  DEFINITION  (OPTIONAL): 

_ □  Etch  using  the  dual  chamber  RIE  reactor _ 

2 1  PATTERN  AlGaN/GaN: 

□  Place  sample  in  sapphire  carrier  attached  with  heat  pump  oil 

□  Place  sample  in  ICP/RIE  chamber  for  5:00  (m:s) 

Cl2  32  seem 

BC13  8  seem 

Ar  5  seem 

RIE  75  W 

ICP  500  W 

Pressure  2mT 

□  Remove  sample  and  clean  with  Acetone  and  Isopropyl 

22  REMOVE  Ni  (OPTIONAL): 

□  Mix  the  following 

HN03  (Nitric  acid)  80  mL 

CH3COOH  (Acetic  Acid)  80  mL 

2  H2S04  (  Sulfuric  Acid)  32  mL 

H20  (AS  DESIRED)  0  mL 

□  1  min  rinse  in  above  solution  or  until  Ni  dissolves 

□  3xDIH20  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

23  REMOVE  SiN  (OPTIONAL): 

□  3  min  BOE:DI(l:l) 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 


149 


Purpose 

Nitride  process 


Mask 

Top  Contact 


Process 

Top  Contact 


Step  # 

Process  Step 

1 

SAMPLE  PREPARATION: 

□  3  min  BOE:DI(l:l) 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

2 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

3 

PMGI  COAT: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

4 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

5 

EXPOSE  1813  with  Contact  MASK: 

□  Align  (  very  critical ) 

□  40  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

6 

1813  DEVELOP: 

□  50  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI  H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

7 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  vellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

8 

DEEP  UV  (  254  nm  ): 

□  200  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 
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Step  # 

Process  Step 

9 

PMGI  DEVELOP: 

□  50  sec  101  Develop  @  500  rpm 

□  40  sec  DI  @  500  rpm 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

10 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

11 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

12 

EVAPORATE  Contact  METAL: 

□  Pre-metal  dip  on  Amonium  Hydroxide:DI  (1:4)  for  1  min 

□  Evaporate  100  A/  1000  A/  500  A/200  A  Ti/Al/Ti/Au  (Ti  on  AlGaN  surface) 

13 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

14 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

15 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 
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Purpose 

Nitride  process 


Mask 

Mesa  Isolation 


Process 

Trench  for 
Contact  Pad 


Step  # 

Process  Step 

CLEAN  SAMPLE: 

□  15  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  15  sec  Isopropyl  Alcohol  at  500  RPM 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Mesa  Isolation  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  50  sec  Develop  with  (1:5)  351:DI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography  :  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 
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Step  # 

Process  Step 

PMGI  DEVELOP: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

DEEP  UV  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI. 

02  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) _ 

PATTERN  AlGaN/GaN  : 

□  Place  sample  in  ICP/RIE  chamber  for  3:15  (m:s) 


ci2 

32  seem 

BC13 

8  seem 

Ar 

5  seem 

RIE 

75  W 

ICP 

500  W 

Pressure 

2mT 

REMOVE  1813: 

□  Heat  1165  to  90°C 

□  15  sec  Acetone  at  500  RPM 

□  15  sec  Isopropyl  alcohol  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

Purpose 

Nitride  process 


Mask 

Contacts 


Process 

Contact 


Step  # 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1818  COAT: 

□  Flood  wafer  with  1818 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1818  with  Top  and  Bottom  Contact  MASK: 

□  Align 

□  70  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1818  DEVELOP: 

□  50  sec  Develop  with  (1:5)  351  :DI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  (  254  nm  ): 

□  200  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 
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Step  # 

Process  Step 

DEEP  UV  (  254  nm  ): 

□  200  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

PRE-METAL  DIP: 

□  1  min  Amonium  hydroxide :DI  (1:4) 

□  4x  DI  full  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

EVAPORATE  Contact  METAL: 

□  Evaporate  Ti/Al/Ti/Au  (100  A/  1000  A/500  A/200  A  Ti  on  surface) 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

ALLOYING: 

□  Set  n2  flow  to  10  mm 

□  Run  RTA  at  900  C  FOR  30  sec  with  no  sample  (dummy  run) 

□  Place  sample  on  RTA  and  purge  with  N2  for  5  min 

□  Run  RTA  at  900  C  FOR  30  sec 

□  Let  RTA  cool  to  -250  C  and  retrieve 

Purpose 

Nitride  process 


Mask 

Top  Contact 
PMGI  VIA 


Process 

Contact  Opening 
and  Reflow 


Step  # 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #3: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Top  Contact  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

EXPOSE  1813  with  PMGI  VIA  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  40  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 
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Step  # 

Process  Step 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  (  254  nm  ): 

□  200  sec  Deep  UV  Exposure  @  1 6mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  60  sec  101  Develop  Immersion 

□  4X  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

DEEP  UV  (  254  nm  ): 

□  200  sec  Deep  UV  Exposure  @  1 6mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  60  sec  101  Develop  Immersion 

□  4X  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

Remove  1813: 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  30  sec  H20  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

REFLOW: 

□  Place  in  250°C  oven 

□  Close  door  and  start  timer 

□  Heat  for  65  sec 

□  Remove  from  oven 

□  Cool  wafer 

INSPECT  WAFER: 

□  Inspect  PMGI. 
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Purpose 

Nitride  process 


Mask 

Bridge 


Process 

Bridge 


Init. 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Bridge  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  40  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

EVAPORATE  Bridge  METAL: 

□  Evaporate  300  A/  4000  A  Ti/Au  (Ti  on  surface,  Au  on  top  of  Ti) 

REMOVAL  OF  BASE  METAL: 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  30  sec  DI  water  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

1818  COAT  &  SLICE: 

□  Flood  wafer  with  1818 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

□  Saw  the  devices 
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C.3  Process  follower  for  Third  Mask  set 


Purpose  Mask  Process 

Active  waveguide  Nitride  etch  (dark  field)  WG  definition 


Step 

Process  Step 

SAMPLE  PREPARATION: 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  1  min  ammonium  hydroxide 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT: 

□  Flood  wafer  with  PMGI  (SF-08) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1805  COAT: 

□  Flood  wafer  with  1805 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1805  with  Contact  MASK 

□  Align 

□  14.0  sec  Exposure  (405  nm)  on  Mask  Aligner  (8.60  settings) 

@  20mW/cm2  optical  intensity 

1805  DEVELOP: 

□  60  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  (  254  nm  ): 

□  95  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP: 

□  50  sec  SAL  101  @  500  rpm 

□  Back  side  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

Step 

Process  Step 

INSPECT  WAFER: 

□  Inspect  PMGI  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  flat  with  a  scrap  wafer  underneath) 

EVAPORATE  Contact  METAL: 

□  Pre-metal  dip  on  Amonium  Hydroxide:DI  (1:4)  for  1  min 

□  Evaporate  250  A  Ti  /900  A  Ni 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

□  Measure  metal  mask  height 

PATTERN  AlGaN/GaN: 

□  Place  sample  in  sapphire  carrier  attached  with  heat  pump  oil 

□  Place  sample  in  ICP/RIE  chamber  for  5:00  (m:s) 

Cl2  32  seem 

BC13  8  seem 

Ar  5  seem 

RIE  75  W 

ICP  500  W 

Pressure  2mT 

□  Remove  sample  and  clean  with  Acetone  and  Isopropyl 

REMOVE  SiN: 

□  3  min  BOE:DI(l:l) 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 
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Purpose 

Nitride  process 


Mask 

Mesa  Isolation 


Process 

Trench  for 
Contact  Pad 


Step  # 

Process  Step 

CLEAN  SAMPLE: 

□  15  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  15  sec  Isopropyl  Alcohol  at  500  RPM 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Mesa  Isolation  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  50  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  #1  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 
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Step  # 

Process  Step 

PMGI  DEVELOP  #1: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

DEEP  UV  #2  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP  #2: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI. 

□  Measure  step  height 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

PATTERN  AlGaN/GaN: 

□  Place  sample  in  ICP/RIE  chamber  for  3:15  (m:s) 

Cl2  32  seem 

BC13  8  seem 

Ar  5  seem 

RIE  75  W 

ICP  500  W 

Pressure  2  mT 

REMOVE  1813: 

□  Heat  1165  to  90°C 

□  15  sec  Acetone  at  500  RPM 

□  15  sec  Isopropyl  alcohol  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 
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Purpose  Mask 

Nitride  process  Bottom  Contacts 

Top  Contacts 


Process 

Contacts 


Step  # 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1818  COAT: 

□  Flood  wafer  with  1818 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1818  with  top  &  Bottom  Contact  MASK: 

□  Align 

□  70  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1818  DEVELOP: 

□  50  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  #1  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

Step  # 

Process  Step 

PMGI  DEVELOP  #1: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

DEEP  UV  #2(  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP  #2: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

PRE-METAL  DIP: 

□  1  min  Ammonium  hydroxide :DI  (1:4) 

□  4x  DI  full  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

EVAPORATE  Contact  METAL: 

□  Evaporate  Ti/Al/Ti/Au  (100  A/  1000  A/500  A/200  A  Ti  on  surface) 

REMOVAL  OF  BASE  METAL: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher 

ALLOYING: 

□  Set  n2  flow  to  1 0  mm 

□  Run  RTA  at  750  C  FOR  30  sec  with  no  sample  (dummy  run) 

□  Place  sample  on  RTA  and  purge  with  N2  for  5  min 

□  Run  RTA  at  750  C  FOR  30  sec 

□  Let  RTA  cool  to  -250  C  and  retrieve 

Purpose 

Nitride  process 


Mask 

Top  Contact 
PMGI  VIA 


Process 

Contact  Opening 
and  Reflow 


Step  # 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

PMGI  COAT  #1: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #2: 

□  Flood  wafer  with  PMGI  (SF-1 1) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

PMGI  COAT  #3: 

□  Flood  wafer  with  PMGI  (SF-11) 

□  30  sec  Spin  at  4000  RPM,  Ramp=200 

□  1  min  Air  Bake 

□  5  min  270  °C  Hot  Plate  Bake 

□  Cool  Wafer 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 10°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Top  Contact  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

EXPOSE  1813  with  PMGI  VIA  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  40  sec  Develop  with  (1:5)  35 EDI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 
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Step  # 

Process  Step 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography:  □  Open  □  Clean  □  Sharp  Definition 

DEEP  UV  #1  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP  #1: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

DEEP  UV  #2  (  254  nm  ): 

□  250  sec  Deep  UV  Exposure  @  16mW/cm2  optical  intensity 

PMGI  DEVELOP  #2: 

□  50  sec  Develop  with  SAL  101  @500  rpm 

□  40  sec  DI  @  500  rpm 

□  Bucket  Rinse  with  DI  H2O 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  WAFER: 

□  Inspect  PMGI.  Under  cut  should  be  obvious  (0.5  pm) 

Remove  1813: 

□  30  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  30  sec  H20  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

O2  PLASMA  ASHER: 

□  4  min  LFE  Barrel  Asher  (place  in  front) 

REFLOW: 

□  Place  in  250°C  oven 

□  Close  door  and  start  timer 

□  Heat  for  70  sec 

□  Remove  from  oven 

□  Cool  wafer 

INSPECT  WAFER: 

□  Inspect  PMGI. 

Step  # 

Process  Step 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

1813  COAT: 

□  Flood  wafer  with  1813 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 1 0°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

EXPOSE  1813  with  Bridge  MASK: 

□  Align 

□  50  sec  Exposure  (405  nm)  on  Mask  Aligner  (90  settings) 

@  20mW/cm2  optical  intensity 

1813  DEVELOP: 

□  40  sec  Develop  with  (1:5)  351:DI  H2O  at  500  RPM 

□  30  sec  DI H2O  Rinse  at  500  RPM 

□  Rinse  wafer  backside  with  DI  H2O  Bucket  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

INSPECT  LITHOGRAPHY: 

□  Place  wafer  flat  towards  top  of  microscope 

□  Inspect  wafer  alignment  with  yellow  filter  on  microscope 

□  Check  Lithography  :  □  Open  □  Clean  □  Sharp  Definition 

EVAPORATE  Bridge  METAL: 

□  Evaporate  300  A/  4000  A  Ti/Au  (Ti  on  surface,  Au  on  top  of  Ti) 

METAL  LIFT-OFF: 

□  Heat  1165  to  90°C 

□  Tape  lift  off 

□  Remove  remaining  metal  using  Acetone  on  a  petri  dish 

□  20  sec  spray  with  Acetone  Bottle  at  500  RPM 

□  30  sec  Isopropyl  Alcohol  at  500  RPM 

□  30  sec  DI  water  at  500  RPM 

□  Dry  with  N2  at  500  RPM 

□  Dry  with  N2  on  clean  Texwipes 

□  Check  backside  and  clean  if  necessary 

DEHYDRATION  BAKE: 

□  1  min  1 10°C  Hot  Plate  Bake 

1818  COAT  &  SLICE: 

□  Flood  wafer  with  1818 

□  30  sec  Spin  at  4,000  RPM,  Ramp=200 

□  75  sec  1 1 0°C  Hot  Plate  Bake 

□  Check  cleanliness  of  wafer  backside 

□  Saw  the  devices 

REMOVE  PMGI: 

□  2  min  1165  at  90°C 

□  6x  DI  H2O  Rinse 

□  Dry  with  N2  on  clean  Texwipes 

APPENDIX  D 


MASK  TEMPLATE  ARRANGEMENT 
D.l  Second  mask  template  arrangement 


D3000-14 

D3030-20 

D3040-20 

D3000-14 

D3030-20 

D3040-20 

D3000-14 

D3030-20 

D3040-20 

D3000-14 

D3030-20 

D3040-20 

D3010-14 

D3020-20 

D3050-20 

D3010-14 

D3020-20 

D3050-20 

D3010-14 

D3020-20 

D3050-20 

D3010-14 

D3020-20 

D3050-20 

D3020-14 

D3010-20 

D3060-20 

D3020-14 

D3010-20 

D3060-20 

D3020-14 

D3010-20 

D3060-20 

D3020-14 

D3010-20 

D3060-20 

D3030-14 

D3000-20 

D3070-20 

D3030-14 

D3000-20 

D3070-20 

D3030-14 

D3000-20 

D3070-20 

D3030-14 

D3000-20 

D3070-20 

D3040-14 

D3200-14 

D3080-20 

D3040-14 

D3200-14 

D3080-20 

D3040-14 

D3200-14 

D3080-20 

D3040-14 

D3200-14 

D3080-20 

D3050-14 

D3190-14 

D3090-20 

D3050-14 

D3190-14 

D3090-20 

D3050-14 

D3190-14 

D3090-20 

D3050-14 

D3190-14 

D3090-20 

D3060-14 

D3180-14 

D3 100-20 

D3060-14 

D3180-14 

D3 100-20 

D3060-14 

D3180-14 

D3 100-20 

D3060-14 

D3180-14 

D3 100-20 

D3070-14 

D3170-14 

D3 110-20 

D3070-14 

D3170-14 

D3 110-20 

D3070-14 

D3170-14 

D3 110-20 

D3070-14 

D3170-14 

D3 110-20 

D3080-14 

D3160-14 

D3 120-20 

D3080-14 

D3160-14 

D3 120-20 

D3080-14 

D3160-14 

D3 120-20 

D3080-14 

D3160-14 

D3 120-20 

D3090-14 

D3150-14 

D3 130-20 

D3090-14 

D3150-14 

D3 130-20 

D3090-14 

D3150-14 

D3 130-20 

D3090-14 

D3150-14 

D3 130-20 

D3100-14 

D3140-14 

D3 140-20 

D3100-14 

D3140-14 

D3 140-20 

D3100-14 

D3140-14 

D3 140-20 

D3100-14 

D3140-14 

D3 140-20 

D3110-14 

D3130-14 

D3 150-20 

D3110-14 

D3130-14 

D3 150-20 

D3110-14 

D3130-14 

D3 150-20 

D3110-14 

D3130-14 

D3 150-20 

D3120-14 

D3120-14 

D3000-14  SO. 5 

D3120-14 

D3120-14 

D3000-14  SO. 5 
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D.2  Third  mask  template  arrangement 


D3000-14 

D3300  14S05R 

D3000  14S10  R 

D3000-14 

D3300  14S05L 

D3000  14S10L 

D3000-14 

D3300  14S05R 

D3000  14S10  R 

D3000-14 

D3300  14S05L 

D3000  14S10L 

D3010-14 

D3300  14S05R 

D3000  14S10  R 

D3010-14 

D3300  14S05L 

D3000  14S10L 

D3010-14 

D3300  14S05R 

D3000  14S10  R 

D3010-14 

D3300  14S05L 

D3000  14S10L 

D3020-14 

D3000  14S05R 

D3300  14S10  R 

D3020-14 

D3000  14S05L 

D3300  14S10L 

D3020-14 

D3000  14S05R 

D3300  14S10  R 

D3020-14 

D3000  14S05L 

D3300  14S10L 

D3030-14 

D3000  14S05R 
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D3000  14S05L 

D3300  14S10L 

D3030-14 

D3000  14S05R 

D3300  14S10  R 

D3030-14 

D3000  14S05L 

D3300  14S10L 

D3040-14 

D3200-14 

D3000  14S15  R 

D3040-14 

D3200-14 

D3000  14S15L 

D3040-14 

D3200-14 

D3000  14S15  R 

D3040-14 

D3200-14 

D3000  14S15L 

D3050-14 

D3190-14 

D3000  14S15  R 

D3050-14 

D3190-14 

D3000  14S15L 

D3050-14 

D3190-14 

D3000  14S15  R 

D3050-14 

D3190-14 

D3000  14S15L 

D3060-14 

D3180-14 

D3300  14S15  R 

D3060-14 

D3180-14 

D3300  14S15L 

D3060-14 

D3180-14 

D3300  14S15  R 

D3060-14 

D3180-14 

D3300  14S15L 

D3070-14 

D3170-14 

D3300  14S15  R 

D3070-14 

D3170-14 

D3300  14S15L 

D3070-14 

D3170-14 

D3300  14S15  R 

D3070-14 

D3170-14 

D3300  14S15L 

D3080-14 

D3160-14 

D3000  14S10  R 

D3080-14 

D3160-14 

D3000  14S10L 

D3080-14 

D3160-14 

D3000  14S10  R 

D3080-14 

D3160-14 

D3000  14S10L 

D3090-14 

D3150-14 

D3000  14S10  R 

D3090-14 

D3150-14 

D3000  14S10L 

D3090-14 

D3150-14 

D3000  14S10  R 

D3090-14 

D3150-14 

D3000  14S10L 

D3100-14 

D3140-14 

D3120  14 

D3100-14 

D3140-14 

D3120  14 

D3100-14 

D3140-14 

D3120  14 

D3100-14 

D3140-14 

D3120  14 

D3110-14 

D3130-14 

D3300  14S10  R 

D31 10-14 

D3130-14 

D3300  14S10L 

D31 10-14 

D3130-14 

D3300  14S10  R 

D3110-14 

D3130-14 

D3300  14S10L 

D3120-14 

D3120-14 

D3300  14S10  R 

D3120-14 

D3120-14 

D3300  14S10L 
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